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INTRODUCTION 


For  geodetic,  ger dynami ca  1  ,  and  astronomical  purposes, 
two  basic  coordinate  systems  a  re  needed:  an  inertial  system 
and  an  earth-fixed  system.  These  two  systems  are  related 
through  precession,  nutation,  and  polar  motion. 

In  order  to  relate  these  two  systems  in  a  precise 
fashion,  we  need  not  only  highly  sophisticated  observation 
methods  such  as  doppler,  lunar  laser,  and  VLB  I  observations, 
but  also  very  accurate  theories,  which  take  into  account  the 
elasticity  of  the  earth's  mantle,  as  well  as  effects  due  to 
the  liquid  core.  A  clear  understanding  of  this  complex  matter 
is  a  prerequisite  also  for  practical  work  in  this  field. 

The  present  report  is  intended  as  a  systematic  review 
which  presents  the  basic  principles  in  a  rather  detailed 
manner  and  is  thus  suitable  as  an  introduction  even  for  sci¬ 
entists  with  little  or  no  previous  knowledge  of  the  field. 

The  report  is  restricted  to  those  theories  which  regard 
the  earth  either  as  a  rigid  body,  such  as  Kinoshita's  recent 
theory,  or  as  a  purely  elastic  solid,  such  as  McClure's  work, 
or  as  a  body  consisting  of  a  rigid  mantle  and  a  liquid  core, 
the  so-called  Poincare  model.  None  of  these  models  is  fully 
realistic,  but  each  contains  important  features  which  form 
a  basis  indispensable  for  understanding,  and  even  for  practi¬ 
cally  and  numerically  treating,  a  more  realistic  model  con¬ 
sisting  of  an  elastic  mantle  and  a  liquid  core.  The  consider¬ 
ations  of  models  of  the  latter  kind,  which  is  a  complicated 
and  difficult  subject, ’will  be  deferred  to  another  report. 
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Emphasis  is  on  the  treatment  of  an  elastic  earth  on  the 
basis  of  Liouvi lie's  equation,  leading  to  a  systematic  tneory 
of  polar  motion,  precession,  and  nutation  of  various  axes 
(rotation  axis,  angular  momentum  axis,  figure  axis,  and  the 
so-called  celestial  reference  pole),  and  on  the  eigenvalue 
problem  for  rotation,  leading  to  a  similar  theory  for  a 
rigid  earth  and  for  the  Poincare  model. 


1.  THE  TIDAL  POTENTIAL 


Precession,  nutation,  forced  polar  motion,  and  earth 
tides  all  have  a  common  cause:  the  gravitational  attraction 
of  sun  and  moon.  Therefore,  the  potential  of  this  attraction 
the  tidal  potential,  plays  a  fundamental  role  in  all  these 
phenomena . 

Consider  the  gravitational  attraction  of  the  moon 
(the  sun  can  be  treated  analogously)  at  a  point  P  on  the 
earth's  surface  which,  to  an  accuracy  sufficient  for  the 
present  purpose,  can  be  represented  by  a  sphere  of  radius  a 
(Fig.  1.1).  The  potential  of  this  attraction  at  P  is 


^  *  Gu  I 


n  =  0  d 


n+  1 


Pn(cosO 


(1-1) 


FIGURE  1.1.  The  tidal  attraction 
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on  expanding  1/1  into  a  series  of  Legendre  polynomials 
P  n  (  c  o  s  ; )  (Heiskanen  and  Moritz,  1967,  p.  33).  The  notations 
are  evident  from  Fig.  1.1;  G  is  the  gravitational  constant 
and  u  is  the  moon's  mass.  We  follow  (Moritz,  1980,  sec.  55). 

The  zero  and  first  degree  terms  (n=0  and  1)  do  not 
cause  genuine  deformations  and  are  therefore  omitted.  The 
dominant  term  is  n=2  ;  h i ghe r- de g ree  terms  are  small  and  will 
be  neglected.  Thus  there  remains 

v  =  v:  =  Gu^r  P.;(cos .  )  .  (1-2) 

This  spherical  harmonic  of  second  degree  will,  in  the  follow¬ 
ing,  be  considered  as  our  tidal  potential. 

Let  us  now  express  P: ( cos v )  in  terms  of  the  geocentric 
spherical  coordinates  of  P  and  of  the  moon's  center.  In  the 
usual  earth-fixed  equatorial  system  the  point  P  has  the  co¬ 
ordinates  where  e  =  9  0  ^  is  the  polar  distance  of 

P  ,  :  denoting  the  geocentric  latitude,  and  >  is  the  geo¬ 
centric  longitude.  Similarly,  the  moon  has  the  coordinates 
(p,h),  where  the  polar  distance  is  given  by  p  =  90  -  i  ,  • 

being  the  declination  of  the  moon,  and  h  denotes  the  Green¬ 
wich  hour  angle  of  the  moon,  that  is,  the  angle  between  the 
Greenwich  meridian  and  the  meridian  passing  through  the 
moon's  center.  Contra  ry  to  astronomical  usage,  both  \  and 
h  are  counted  positively  towards  east  (Fig.  1.2). 

Then  P  ■  ( c o s  v  )  can  be  expressed  by  means  of  the  decom¬ 
position  formula  for  spnerical  harmonics  (Heiskanen  and 
Moritz,  1967,  p.  33),  and  we  obtain 


V  2  -  V  2  o  +  V  2  ;  +  V22  > 

with  the  zonal  part 

v 2 o  =  Gu  Jr  P2g (cose )P20 (cosp)  , 

the  tessera!  part 

i  : 

v  2 }  =  -j  Gm  P;i(cose)P2i(cosp)  cos('-h) 

and  the  sectorial  part 

v  -  2  =  -pr  G  u  Pj?  (cosh  )  P  2  ->  (cosp)  C  0  s  2  (  •.  -  h  ) 


Tne  Legendre  functions  are  given  by  the  well-known  expressions: 

P  •  ( cos  -  ) 

P,  : (cos-) 

The  coordinates  of  the  moon,  p  and  n,  are  functions 
of  the  time  since  the  moon  moves  along  its  orbit.  Tnerefore  , 
eqs.  (1-4)  to  (1-6)  can  be  considered  spherical  harmonics  in 
and  ■  with  coe f f i r  \  e n ts  that  are  functions  of  time.  We 
expand  these  functions  into  trigonometric  series  and  write 
the  result  in  the  form  (Doodson,  1922;  McClure,  1972,  p.  92): 


P  ( c  O  S  •  )  =  ?  cos 


3  s  i  n  c  o  s 


P  (cos-)  =  3s  in- 


v  =  v 


v 


nm 


+  V  ’  +  V 


(cos-) 


(  1  -3 ) 


•  A  .cos  !  j  .  t  +  .:  ..  +  m  '•  +  (n-m)*‘ 

n njd  [_  n m  j  u  n m  ]  i  '  Z 

j 

Here  n  =  2;  m  =  0 , 1  ,  2  ;  the  index  d  numbers  moon  ( d  =  1 )  and  sun 
(d-2)  ,  ...  denotes  the  masses  of  moon  and  sun,  and  c  de- 

notes  tne  mean  radii  of  the  lunar  and  the  solar  orbit  (con¬ 
sidered  with  respect  to  the  earth).  As  an  example,  we  write 
( 1  -  3 )  explicitly  for  the  case  m  = 1  ,  which  will  be  of  parti- 
cular  importance: 


G..  . 

— 1  a  P.  : ( cos  -  )  • 

°  i  (  1  -  9  ) 

A  ■.  •  s  i  ri  (  j  -  ■  ,  j  t  + 


+  '■  ' 


since  c  o  s  (  :i  +  -  /  2  )  =  - s i n  i  .  Omitting  the  subscripts 
in  the  coefficients  and  considering  the  effect  only  of 
the  moon  (or  only  of  the  sun)  we  have 


'2  1  = 


G  u 


c 


a  P  -  ,  (  c  o  s 


)  A  s  i  n  (  „  .  .  t 


(1-10) 


This  simpler  expression  will  frequently  be  used  later. 1 

We  remark  that  the  arguments  are  linear  combinations 
(with  integer  coefficients)  of 


s 

h 

P 

N 


that  is 


lunar  mean  longitude, 

solar  mean  longitude  (not  to  be  confused 
with  the  hour  angle  as  used  above), 
mean  longitude  of  lunar  perigee  (same  re- 
ma  r  k  )  , 

longitude  of  the  mean  ascending  node  of  the 
lunar  orbit, 

mean  longitude  of  solar  perigee, 
local  mean  lunar  hour  angle. 


>  .  t  +  s  ,  +  m\=n,T  +  n-s  +  n,h 

nm] a  nm] d 


nup  +  n  r  N  +  n.Ps 


(  1-U) 


(Melchior,  1978,  p .  33;  McClure,  1973,  p.  95). 


l)  Mote  that  the  sum 


is  an  infinite  s e r 1 e s ! 


Well  known  are  the  tidal  de vel opements  of  Doodson  (1922 
and  Cartwright  and  Tayler  (1971);  see  also  (Cartwright  and 
Edden,  1973).  A  new  developement  has  been  given  by  Heikkinen 
(  19  78  )  . 


2.  ROTATION  OF  A  RIGID  BODY 


The  basic  equation  for  the  rotation  of  a  rigid  body 
is  very  simple: 


the  derivative  of  the  angular  momentum  H  with  respect  to 
time  t  equals  the  torque  L  ;  both  H  and  L_  are  vectors, 
which  is  indicated  by  underlining  their  symbols. 

This  equation  holds  in  a  nonrotating  (inertial)  coordi¬ 
nate  system;  in  a  body-fixed  coordinate  system  (which  rotates 
with  the  body),  its  equivalent  is 

3H 

3t  +  Id.  H  =  k  •  (2-2  1 

Here  3 / 3 1  denotes  the  time  derivative  in  the  body-fixed 
system  and  w  is  the  rotation  vector  whose  direction  coinci¬ 
des  with  the  instantaneous  axis  of  rotation  and  whose  magni¬ 
tude  is  the  angular  velocity  w  of  rotation;  the  cross  (■) 
denotes  the  vector  product  of  two  vectors.  These  equations 
can  be  found  in  any  text  on  mechanics;  cf.  (Synge,  1960). 

Eq  .  (2-1)  is  fundamental  for  precession  and  nutation, 

which  is  the  motion  of  the  earth's  axis  in  inertial  space; 
and  (2-2)  is  basic  for  polar  motion,  which  is  the  motion 
of  the  earth's  axis  with  respect  to  an  earth-fixed  coordinate 
system. 
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The  velocity  v  of  a  point  of  tne  body  is 

v  =  ■  x  (  2  -  3  ) 

x  denoting  the  position  vector  of  the  point.  This  relation 
is  substituted  into  the  equation  defining  the  angular  momen¬ 
tum  H  , 

H  =  ///x  «  vdM  ,  (2-4) 

in  which  the  integration  is  over  the  body,  dM  being  the 
mass  element.  The  result  is 


H  =  C 


(2-5) 


where  C  is  a  tensor  (a  3x3  symmetric  matrix 


inertia  tensor.  The  elements 
gi  ven  by  the  formul a 


C 

3-D 


of  the  matrix 


,  the 
C  are 


C  .  . 
i  D 


///(*, 


k  '  i  j 


x  .  x  ) dM 


(2-6) 


using  index  notation:  i  ,  j  ,  k  run  from  1  to  3  ,  x  •,  =  x, 
x  =  y,  X;  =  z  are  coordinates  in  a  body-fixed  Cartesian 
System,  5  =  1  if  i=j  and  0  if  i/j  (Kronecker  delta 

and  summation  over  repeated  subscripts  in  a  product  (:<  in 
the  formula)  is  prescribed;  cf.  (Jeffreys,  1931). 

Equivalent  but  somewhat  more  explicit  is  the  form 
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J 

x 

-D 

A 

-D 

X 

c  - 

-D 

x  y 

J. 

'0 

-D 

-D, 

•- 

J  _ 

where  the 

diagonal 

elements 

a  re 

mo  me  n  t  r> 

o  t  i  ne  r  t  ’  a  ,  n 

J 

X 

=  ;7/'y-  + 

?  • 

)  dM 

’ 

and  the  off-diagona 

1  terms 

are 

products 

of  inertia,  e 

D 

XV 

=  / J/xydM 

If 

the  principal 

axes 

0  f 

i  ne  r  1 1  a 

are  chosen  as 

nate  axes 

,  then  the 

tensor 

C 

assumes 

diagonal  form: 

A 

0 

0 

|n 

ii 

o 

B 

o  ;  . 

;  0 

0 

C  i 

A ,8 ,C  being 

the 

p  r  i  n  c  i  p  a  1 

moments  of  inertia. 

In  this 

case ,  (2-5) 

reduces  to 

q  '  v  e  s  l  'j  I  e  >'  s  equations 


f.  l  i  v  >'  t  a  *  •  •  j  /  r  ■  .  p  i  mo  turn;  t  n  e  time  a  e  r  i  v  a 
1  ►-  n  t  ■  •  >■  ■;  /  i  t  t  . 


'  .1 1  1 

>  n  j  r  i 

1  •  i  e  i  r  r  n  .  e  t  u  : 

-  cive  a 

v  £  r'  / 

•1  f  1 . 

n  n  :  *  m  t *  •  t?  * 

■  1  n  >,  t  v  ;;0  1  i  •'  mo  1 1  o  r, 

}  *  *  n  £ 

’ll-  •; 

r,  . .  r  .  i 

l.  e  x  e  o  t  e  'J 

"j  4  tun  a 

no  ^ c  :  r 

r  a  t  i 

n  a  !  ;  /  m  in  •*  r  i 

■  /  v  ’’Ut  B  -  A  . 

'nen 

1  ^ 

J.l  '  i  ■  '  r.  I've 


mine  <1  1  .1  ?  e  / 


-Just. 
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we  transform  the  first  and  second  equation  into 


of  which  the  solution  is 

u;  =  iCOS(u  t  +  ,  )  , 

W2  =  *  S  in(  ;£t  +  ,  )  ,  2 " 1 0  ; 

with  constants  i  and  ,  .  From  (2-16;  we  have 

ij]2  +  =  const.  ,  '2-17) 

which  is  the  equation  of  a  circle.  Together  with  *  ,4  =  c  o  n  s  t . 
this  means  that  the  rotation  a<is  describes  a  circular  cone 
around  the  axis  of  symmetry  (Fig.  2.1  .  The  angle  of  aper¬ 
ture  is  about  0.2";  the  period  T  is  obtained  from 
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(2-14)  as 


T  --  ~  --  --  -  305  days  .  (2-13) 

L  -  A 

E 

T  h  i  s  E  uj  e  r  p  e  r  ijod  would  hold  if  the  earth  was  a  rigid  body. 

Tne  fact  that  toe  actual  period,  tne  Chandler  period ,  i s  about 
430  da  vs,  indicates  that  the  earth  is  not  rigid. 

Regarding  terminology,  tne  constant  -  defined  o y 
(2-14'  will  be  called  the  Euler'  frequency ;  it  will  play  a 
basic  role  throughout  the  present  report. 

Strictly  speaking,  is  an  "angular  frequency",  where¬ 

as  the  name  "frequency"  is  usually  reserved  for  the  quantity 


?  - 


1 


now ever,  we  shall  consistently  speak  of  frequency  in  the  sense 
of  "angular  frequency". 

The  coefficient  i  in  (2-16)  is  called  amp  1 i tude  , 
ind  .  is  the  phase. 


1  3 


3 .  THE  LIOUVILLE  EQUATION 


The  basic  equation  (2-1)  , 
dH 

dT  =  —  *  t3'1' 

is  really  quite  general.  It  holds  for  the  rotation  of  an  ar¬ 
bitrary  body,  rigid  or  not  (Truesdall  and  Toupin,  1960,  p .  5  3 1  )  . 
T'he  underlying  coordinate  system  is  an  inertial  system  which  we 
denote  by  XYZ  or  X  ;  X  :  X  -  . 

In  a  rotating  system  xyz  or  x  ;  x  x  .  which  is  attached 
to  the  rotating  nonrigid  body  in  a  way  to  be  explained  later, 
the  angular  momentum  equation  takes  again  the  form  (2-2), 


3H 

—  +  'j  x  H  =  L  , 
3  t 


but  ( 2-5  )  is  genera  1 i zed  : 

H  =  C  w  +  h  ( 3-3  ! 

where  n_  is  a  relative  angular  momentum  defined  by 

h  =  ///  x  *  u  dm  .  (  3- 4 ) 

Here  u  is  the  velocity  with  respect  to  the  system  x;x:x:, 
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which  is  related  to  the  velocity  v  with  respect  to  the 
inertial  system  X:X_.X-,  by 

v  =  ^  ■  x  +  u  .  (3-5) 

Equations  (3-3)  and  (3-4)  are  readily  obtained  by 
substituting  (3-5)  into  the  defining  equation  (2-4). 

The  inertia  tensor  C  is  again  given  by  (2-6). 

The  meaning  of  these  equations  is  easily  understood. 

If  the  earth  (our  rotating  body  will  always  be  the  earth) 
is  not  rigid,  then  there  is  no  coordinate  system  at  which 
all  particles,  of  which  the  earth  is  composed,  are  at  rest. 
Thus  they  move  with  respect  to  our  system  x;  x;x3  with 
velocity  u_  ,  which  is  considered  small  since  it  is  zero  for 
a  rigid  body.  Thus  (3-5)  differs  from  (2-3)  by  a  non¬ 
zero  u.  This  relative  velocity  u_  causes  the  relative 
angular  momentum  (3-4)  to  be,  in  general,  different  from 
k  zero. 

By  substituting  (3-3)  into  (3-2)  we  get 

^■(C1+h)+y<(C_o><h)  =  L,  (3-6) 

which  is  called  the  L i ouvi 1 1 e  equation  (Munk  and  Macdonald, 
1960,  p.10).  This  equation  will  be  fundamental  for  the  mathe 
matical  description  of  polar  motion  for  a  nonrigid  earth. 


5 


It  is  now  of  basic  importance  that  the  axes  x  x  x 
can  be  chosen  such  that  h  =  0 .  They  have  the  property  that 

///  u  :dM  =  mi n 1  mum  (3-7  } 

(Jeffreys,  1970,  sec.  7.04)  and  are  called  Tisserand  axes 
(Munk  and  Macdonald,  1960,  p.10).  Then  the  basic  equations 
are  formally  the  same  as  for  a  rigid  body,  eqs  .  (2-2)  and 
(2-5),  but  note  that  now  the  inertia  tensor  C,  eq .  (2-6), 
will  be  a  function  of  time  since  the  shape  of  the  body  will, 
in  general,  change  with  time.  In  the  following  we  shall  al¬ 
ways  use  Tisserand  axes. 

This  is  convenient  as  long  as  one  disregards  relative 
motions  such  as  ocean  currents  and  winds,  as  we  shall  do. 

For  the  consideration  of  such  effects  see  (Munk  and  Macdonald, 
1960,  p.  123;  lambeck  and  Cazenave,  1973,  1974;  Capitaine, 
1980;  Lambeck,  1980). 

Linearization.  For  Tisserand  axes,  the  Liouville 
equation  (3-6)  may  be  written 

yt  )  +  -  x  w  )  =  L  •  (3-8) 

We  shall  now  linearize  this  equation  as  follows  (Munk  and 
Macdonald,  1960,  p.37;  McClure,  1973).  The  inertia  tensor 
is  wr i tten 

C  =  C0  +  £  (3-9) 


i 
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where 

[ 

!  A 

0  0  ! 

i 

C  ;  ; 

c  .  .  C  ;  : 

C,  - 

!  o 

A  0  :  , 

i 

c  =  : 

c 

c.  .  c:?  .  (3-10) 

1  0 

0  C 

i 

c  ;  ? 

C  _  .  C  :  ; 

Thus, 

C  3 

corresponds  to 

the  model 

o  f 

an  undeformed  earth 

whose 

principal 

axes  of  inertia  coincide  with  the  coordinate 

axes  , 

which 

has 

ro  ta  t i ona  1 

symme  t  ry 

(E 

-A)  and  whose  princi- 

pal  moments 

0  f  1 

n  e  r  t  i  a  A 

and  C 

a  re 

cons  tant  in  ti me  .  The 

tensor 

c 

takes 

into  account  the  deviation  of  the  actual 

earth 

from 

this 

s  i  in  p  1  i  f  i  e  d 

model  . 

The  rotation  vector  ^  is  written  as 


£  +  6,  (3-11) 

where 

'■  0  ' 

=  |  0  ,  (3-12) 

corresponds  to  a  rotation  with  constant  angular  velocity 
around  the  z  -  a  x  i  s  and 


m, 

in  .  (3-13) 

m  . 


expresses  deviations  of  the  rotation  axis  from  the  z-axis 
( m  ■  and  rn  _  )  and  variations  of  the  rotational  speed  :  :r, 
Both  c  and  are  considered  small  quantities  whose 

squares,  products  and  higher  powers  can  be  neglected. 

We  substitute  (3-9)  and  (3-11)  into  r  3-3 ;  and 
retain  linear  terms  only.  The  result  is 


L  . 

=  A 

;.!ti .  +  (  C  -  A  ) in  +  :.t  -  =  - 

■  ■•'c.  .  , 

L; 

=  A 

,i(ti  -  (  C  -  A  )  -  m  ;  +  ...  0  ;  ;  + 

- C  .  .  , 

(3 

L  ; 

=  c 

Itl  :  +  t 

the 

dot 

denoting  a  time  derivative 

•  /  i  t . 

1 1 

will  be  very  convenient  to 

combine  the  first 

two 

eq  ua 

t  i  o  n  s  by 

using  complex  notation. 

D  u  1 1  i  h  q 

m 

=  m  ■_ 

+  i  m:  , 

c 

=  C  ; 

;  +  i  C  ;  3  , 

■  m 

L 

=  L; 

i  L 

where  i 

-  = 

-1  (note  the  difference 

between 

the 

comp  1 

e  x 

be  r 

L 

and 

the  three-vector  L ,  and  between 

the 

comp 

1  e  x 

number 

c  and  the  tensor  c  i).  The 

result  i 

s 

L 

=  A;. 

it  -  i  (  C  -  A  ) m  +  :.t  +  i ..  • 

c  ; 

(  3 

i  n  f 

act. 

by 

substituting  (3- IB),  performing 

the 

comp 

I  e  x 

mult 

i  p  1  ’ 

cations  and  separating  real 

and  i ma g i 

na  ry 

part 

s  , 

qe  t 

back 

the 

first  two  equations  of 

(3-14). 

I 


E q .  (3-16)  can  be  written  in  t n e  form 

rii  =  i  • ..  (  m  -  .  } 

w  here 

C-_A 

E  "a 


•s  again  the  Euler  frequency  (2-14;  :we  may  identify  tne 

constant  ...  in  sec.  2  with  the  present  )  ;  and 

_ U. _ +  _c_  _  1 1 

(  C  -  A  }  ;.  -  C-A  (  C  -  A  )  :,  ;3'19 

Tne  o u  a  n  t i t  y  -  is  called  tne  excitation  function  because 
it  causes  a  deviation  of  polar  motion  from  the  simple  case 
of  free  rotation  of  a  rigid  earth  considered  in  sec.  2,  for 
w n i c  h  . =0 . 

If  the  function  .  is  known,  then  (3-17)  can  be  sol¬ 
ved  for  m ,  obtaining  tne  components  m  and  m  character 
icing  a  deviation  of  the  earth’s  rotation  axis  from  tne  c  -  a  x  i 
that  is,  polar  motion.  Similarly,  the  third  equation  of  (3-1 
w n  i  c h  can  be  written  a s 


can  be  solved  for  m  to  get  variations  of 
f  ation,  )  r  v  .i  r  i  >.  t  i  on  s  of  t  h  e  1  e  n  a  t  n  of  day. 


the  speed  of  r o - 


In  other  terms,  polar  motion  i;  the  variation  or  tne 
direction  of  the  rotation  vector  and  m  c h a  r a c  t e r i  ; -  ; 

variations  in  the  l_e  nqt  h  of  .  .  It  is  ver  /  r  e  mars  a  u  !  *.  ■  <  t 

both  phenomena  are  separated  in  the  linear  approx i mu  1 1 •  n . 
the  sequel  we  shall  restrict  ourselves  to  polar  mu  t  ion  ,v  n  •. 
is  more  interesting  and  less  simple. 

We  shall  also  linearize  the  expression  f  3  -  3  ;  ,  wn:cr 
for  h  =  0  reduces  to 

H  =  C  .  '3-2  1 

The  substitution  of  (3-9)  and-  (3-11)  yields,  on  neglect! 
terms  of  second  and  higher  order, 

H  =  A m  f.  c  ,  '3-32 

where 

H  =  H-  -t-  iH^  ,  (3-22 

and 

H  3  =  C +  C.tm-j  +  f.  Ci  :  .  (  3  -  2  ~ 

The  division  of  the  equatorial  components  H;  and  H-  of 
the  vector  H  by  its  length,  which  approximately  is  C., 
gives  the  equatorial  components  h,  and  h  of  the  unit 


2  0 


vector  H H  : 

,  H  ■  ,  H .  u 

h  ;  -  p-- ■  ,  h  ”  c~"  ’  n  =  n  •  t  1  h ;  .  [  , 

These  quantities  define  the  deviation  of  the  angular  mo men - 
turn  axis  (the  direction  of  hi)  from  tne  z-axis,  in  the  same 
way  as  in  and  m  _  define  tne  d  e  v  1  a  t  i  o  n  of  the  rotation 
axis  from  the  z-axis.1 

Finally,  we  consider  the  f i gu re  axis  which  is  tne  axis 
of  maximum  inertia  for  the  deformed  earth,  that  is,  a  prin¬ 
cipal  axis  of  cne  tensor  C  (the  z-axis  as  a  principal 
axis  of  the  tensor  C .  The  determination  of  the  principal 
axes  of  a  symmetric  matrix  is  s t ra i gh t f o rwa rd ;  for  tne  tensor 
! 3 - 9 )  we  get 


C-A  C-A 

or  in  complex  notation 

f  =  f;  +  if.:  =  £T7f  (3-27) 

(McClure,  1973,  Appendix  E).  The  quantities  f  ;  and  f : 
are  equatorial  components  o f  the  unit  vector  of  the  figure 

axis. _ _ _ 

H  Since  the  relative  angular  mo men  turn  vector  (3-4)  is  zero 
the  letter  h  is  free  for  our  present  new  use. 


4.  FREE  POLAR  MOTION  FOR  AN  ELASTIC  EARTH 


If  the  earth  rotates  about  an  axis  which  deviates 
from  the  axis  of  symmetry,  there  occur  centrifugal  forces 
which  tend  to  distort  it,  and  an  elastic  earth  yields  to 
this  distortion  (this  is  similar  to  distorting  forces 
acting  on  an  unbalanced  wheel).  It  is  well  known  that  this 
distortion  introduces  products  of  inertia  c-:  and  c  _  . 
which  are  proportional  to  the  deviation  of  the  rotation 
axis,  m;  and  m?  (Jeffreys,  1970,  sec.  7.04;  Munk  and 
Macdonald,  1960,  p.  38).  In  complex  notation,  using  (3-15), 
we  have 

c  =  £  (C-A)m  ,  (4-1) 

s 

where  k  is  a  Love  number  well  known  from  elasticity  theory, 
and  k  is  an  abbreviation 

s 

k  =  — tlC-A-).  (4-2) 

s  a 5  ■' 2 

and  bears  the  somewhat  unfortunate  name  of  secular  Love 
number  (Munk  and  Macdonald,  1960,  p.  26).  Here  G  denotes 
the  gravitational  constant  and  a  is  the  earth's  equatorial 
radius  or,  to  the  same  accuracy,  its  mean  radius;  a  =  R  = 

=  6371  km.  More  about  the  constants  k  and  \  will  be  said 
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at  the  end  of  this  section. 

This  quantity  (4-1)  is  to  be  used  in  the  excitation 
function  (3-19).  In  the  case  of  free  rotation  to  be  con¬ 
sidered  in  this  section,  there  is  no  effect  of  sun  and  moon, 
that  is,  no  external  torque  L.  Thus  (3-19)  reduces  to 


c  id  k  , 

*  =  "  JUTJq  =  k  {m 

'  s 


(4-3 


the  subscript  "rd"  denoting  "rotational  deformation",  and 
(3-17)  becomes 


rfi  =  i  r  m 
c 


(elastic  earth)  . 


(  4  -  4  ; 


Here 


'  c 


E 


1 


(  4  -  5  ) 


is  the  Chandler  f requency  ,  the  name  will  be  explained  below; 
note  that  n  is  the  Euler  frequency  (2-14). 

The  complex  equation  (4-4)  splits  up  into  two  real 

ones  : 


rf)  ■  +  ;  {  m  =  0  , 

rti-  -  >:  m,  =  0  , 

c  ■ 


(4-6  ) 


which  have  the  same  form  as  (2-15)  which,  in  our  p  re  s  <_•  n  t 
notation,  could  be  written  as 

rt)  =  i  c  „  tn  (rigid  eirtn)  .  (4-7 

With  the  numerical  values 

k  4  0.30,  ks  =  0.96  (dimensionless)  '4-3 

we  ge  t  from  ( 4  -  5  ) 


which  means  that  the  period  (2-18)  is  lengthened  by  the 
factor  1/0.7  4  1.4  .  The  multiplication  by  this  factor  brings 
the  Euler  period  of  305  days  (which  would  hold  for  a  rigid 
earth  and  corresponds  to  the  Euler  frequency  -  )  close  to  tne 
actual  Chandler  period  of  about  430  days,  which  corresponds 
to  the  Chandler  frequency  ;  . 

The  solution  of  (4-4)  may  be  written  in  c o m p 1 e a  form 
as 

m  =  m 0 e  '  : c c  (4-10) 

with  an  arbitrary  (complex)  constant  m?.  This  is  immediate¬ 
ly  verified  by  substitution  into  (4-4),  which  again  shows 
the  advantage  of  complex  notation.  Of  course,  (4-10)  is 
equivalent  to  (2-16),  with  i  replaced  by  -  ( „  . 


A  t-  Ill  ,i  r‘  r  OH  L  J  Vf  r.  u  Hi  l)  lj  f  .  r,  h  .  i  j  .  ^  'j  V  *ir  '  u  iT. 

.  jcdO1-’  ’  !» e  *•  i  a  s  *  i  re  •;  p  o  •,  s  e  o  f  \  '■  e  e  a  r  t n  '  n 

[nj'.  a  r:  e  *  r .?  ••  n  i  c;  l  •>  ti.fl.  i  n  9  potent  1  a  !  v  wrier,  •  •:  3  uppr  e 

to  be  a  s  p  n<-  >•  1 ..  j  !  n  a  rmun  ;  c  'am.  t  ion  of  record  decree,  cdube 

jn  e  1  a  -.  t  1  c  O’  ;  tort  ion  0 1  tne  ea  rtn  and  thereto  a  t'o'tr 
1  t ;  g r a  /  i  t a  t  ;  ;•  n a  ;>u*entia  I  t, / 


it  1  :  J  ■  <,  ;  called  a  1 1  0  A  ;  -e  ff«?  c  t  ;  v  e  _  c  v  e  n  umDej*  V  r  3  r  c 

Ma  0  do  n  3 7  4  ,  iso  0  ,  p  .  /  ,•  . 

7  0  e  U'coiar  Love  number  ►  as  Ge  fined  by  -1-2  , 

on  tne  0  the  r  nano,  "can  be  1  n  te  rp  re  t  ed  as  a  measure  z f  t”-e 

e  a  t  ’■  1  ■,  i  e  1  d  to  centrifugal  deformati  on  in  tne  cause  0  f  it 
<1  eve  i  0  potent  .luring  the  last  five  billion  years  or  so"  •  2_b  1 
p .  25  .  Tnis  is  a  permanent  deformation  wnicn  is  typically 
none  1  a , 1 1  .  :  it  corresponds  to  tne  rotational  deformation  of 
a  f 1 u 1 d  eartn,  peing  numerical]/  equal  to  tne  "fluid  Love 
number "  k ,  ;  1  bid.  ,  p  .  26). 

'ne  difference  between  *■  and  k_  according  to  d - 
may  be  explained  by  a  number  of  h  y  d  0 1  h  e  s  e  s  .  but  it  is  not 
known  wrich  of  then  e  'if  any)  is  c  0 r ret t  .  i  bid.  ,  p  .  2  ~  1  .  ”  n 
usual  approacn  is  to  use  the  elastic  Love  number  k  for 
all  elastic:  de  formations  ,  whether  they  are  caused  by  centri 
fuqal  or-  tidal  disturbances  and  whether  tney  represent  con¬ 
stant  or  temporally  variable  deformations  (cf.  Melchior,  1? 
Recently,  however,  some  authors,  following  McCljre  11973), 
are  using  the  usual  Love  number  k  for  temporally  variable 


de  f  0  rma  t  ’  an  $  a  n  -i  tne  s  e  c  u  1  w  r  L.jvh  numbe  <•  ►  f  o  >■  .  un  ,  t  a  n  \ 
elastic  d  e  r  o  r  m  a  1 1  o  n  s  . 

In  tne  ODiniun  of  the  present  dutnj r ,  'ucn  an  jpproj.n 
is  inaopeopnate.  The  theory  of  de  forna  t  i  on  •,  of  an  elastic 
earth  makes  no  difference  be  tween  deformations  which  are 
constant  in  time  and  for  those  which  ire  temporally  variable 
for  ootn,  the  Love  n  u  in  d  e  r  k  i  s  r  e  1  e  v  o  r  t  .  Even  for  earth 
models  with  a  liquid  core,  fur  which  is  a  slight  dependence 
of  k  on  frequency,  give  for  constant  deformations  a  k 
that  is  close  to  0.3;  cf.  (Wahr,  19/9;.  It  would  require 
an  earth  model  of  a  completely  different  rheoluqy  to  get  for 
constant  deformations  a  k  close  to  1,  if  this  is  at  all 
possible  and  physically  meaningful.  No  such  model  has  oeen 
given  so  far. 

Since  geodynamical  computations  should  De  based  on  a 
well-defined  meaningful  model,  it  is  strongly  advocated  to 
use  a  k  following  from  such  a  model  (Molodensky,  1961; 
Wahr,  1979).  No  physically  observable  error  ’  s  introduced  in 
this  way:  even  if  the  secular  Love  number  were  "true"  for 
constant  deformations,  this  would  only  imply  a  reference 
model  of  a  slightly  different  flattening  but  not  change  at 
all  the  observable  physical  situation. 

Therefore  we  shall  use  k  to  characterize  the  elastic 
response  of  the  earth  also  for  the  constant  part  of  the  de¬ 
formation,  restricting  kg  purely  to  its  use  as  the  abbre¬ 
viation  (4-2)  without  attempting  a  physical  i n terp re ta t i on 
of  k  ^  in  the  sense  of  an  elastic  or  nonelastic  response. 


We  finally  note  tnat,  in  a  rigid  oody,  tnere  are  no 
elastic  deformations  and  the  shape  of  tne  body  and  its  gra¬ 
vitational  potential  V  do  not  change.  Hence  • V  =  0 ,  so 
that  (4-11'  implies 

k  =  0  for  a  rigid  body.  (4-12 

Therefore,  formulas  for  a  rigid  earth  can  be  ODtained  from 
those  for  an  elastic  earth  by  simply  putting  k  =  0 . 


INFLUENCE  )r  SUN  AND  MOON  ON  P';i  *1 0  *  ItN 


he  o a s  i  c  equations  are 


ih  =  i  • .. ;  m  -  .  )  , 


i  n  -j  ;  -  r 


C  1  L 

c^a  "  rc-'Av 


'ne  excitation  function  may  be  s  o  I  •  t  o  as  *■;!!.* 


of  which  the  terms  represent: 


rota  ti  ona  1  d e  f  o  rma  1 1  ) n 
distortions  . 


...  e  f  f  e  c  t  o  r  t  ne  tons- 
moon  , 


e  <  --  •  e  c 


tidal  de  *  o  rma  *  •  n  o  ‘  t  "  -•  -  ir*  n 


"he  ro  t  a  t  i  on  a  1  d  e  c  c  r  n’  a  f  i  a  n  has  area-’,  :<e  -  v. 
preceding  sect  •  an,  .  >hm  t  :  '/on  r  ,  d  - 


two  terms  will  be  treated  now. 


n  i  o  1  a  r  torque.  The  torque  exertec  Dy  the  moon 
irtn  is  equal,  hut  with  opposite  sign,  to  the 
\e-‘teo  by  tne  earth  on  tne  noon  (Kaula,  1963,  sec. 

•;  h  i  )  r ,  19  73,  p.  39).  The  la  ter  is  equal  to  x  •  K , 
;  tne  position  vector  of  the  moon  'Fig.  3.1, 

,  t  ••  tore-;  or  attraction  of  the  moon  oy  tne  earth 
»-  f  i  n  •  t  i  on  o  f  the  vector  product,  and  in  view  of  tne 
a  1  .  ,  'pp  t  rv  of  tne  earth,  tnis  torque  n  as  tre  ma  ;r  ■ 

,v”.  fe  ■  is  tt.e  mer’d'.  an  com;,  oner  »  q  *  •'  ; 


moon 


<_•  I 


2  9 


the  earth,  which  for  rotational  symmetry  tires  tne 
known  f o  rm 


GM 

r 


J 

r 


D  / 


!  C  0  ' 


Of.  !  H  e  i  skanen  and  Moritz,  1967.  o  .  7  3 

ves  to  second -order  spherical  narwonii. 

’  s  t r  e  invitational  Lon/j"t  and  M 
a  -  d 
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where  the  subscript  d  i  s  t  i  n  g  u  i  s  n  e  >  >'jh  and 

finally  the  simpler  complex  fo'-m 

L  =  -  3Gc  -  "  A  e":  ; 

i  -  c  .  ■  i 

as  usual,  L  =  L-.  +  i  L  _  . 

Now  (  5  -  2  )  gives  .  .  ,  using  5  -  r  i  : 


i  L 

TTrr~: 


=  i  •  0  e 


where  the  dimensionless  coefficients  3  are  define: 


B  =  -  3  — V  A  . 

1  .. "  C  ;  i 


Tidal  de  fo ma  t  i  on  .  The  products  of  inertia, 
and  C •  .  ,  are  well  known  to  be  re  1 atea  to  the  tes.er 
In  =  2,  n=l)  o*  the  spheri cal -harmoni c  expansion  c  *  toe 
t  a  t ' o  n  a  1  potential  V  ;  c  f  .  ( H  e i s  k a  n  e  n  and  M  o  r  i  t  a  .  1  1  b 
In  the  present  case,  where  c  •  :  and  c  represen  t  t 
t  j  rbances  o*  the  inertia  tensor  due  to  the  tidal  c  e  *'  o  •• 
tn e  relevant  potential  is 


w  n  '  c  ■!  i  i  tne  c  n  a  n  a  e  of  t  .a  e  n  i  /  .  :  a  t  o  oa  ■  p  o  t  e  r,  a  :  a 

to  the  de  to  :-i  3 1  i  ■  i  n  of  the  earth  ty  tne  t'.  da’  patent',  a’, 

c  f  .  4-11;.  in  v  we  thus  need  only  tne  tessera!  par 

,1-10!,  and  ii'jinn  ••5-14  ,  we  get  in  a  ratner  s  t  r  a  '  g  n  t 
w  a  r  a  w a  y  !  McC  lure,  19  7  3  ,  p  .  22): 

-  r  ..  n  ■;  i  n  (  .  _  t+  .•  . )  , 

-  •  -  '1  cosL  t  +  /  i  , 

:  L  1  "  1  :  ’ 

wnicn  anain  -.an  be  simplified  Py  complex  notation 

c  =  -  i  .  0  -  A  »  B  e~  i  '  • 

r.  i 

here  we  nave  used  ( 4  -  2  j . 

Tidal  deformation  causes  a  change  of  products  of 
e  ••  t :  a  ,  so  5 -  2  )  reduces  to 

c  i  c 

"  =  C-£  ‘  TT-lTf '  ' 

sin,.. (■•  I.  nos  Peon  considered  in  (5-13).  The  aifreren 
•'i  f  '5-1  7  ;  <3  i  v  e  s  c  ,  which  togetner  w  i  t  r  (5-17)  is 
to  ted  into  5  -  l  h  )  .  The  result  can  be  written  as 


where 


'5-20 


and  s.  is  defined  by  (-1-2);  B  is  <3 :  v  e  n  by  «  5  -  1 4  >  . 

3  _1 

The  physical  meaning  of  ..  and  .  snould  be 

1  j  i  l.,‘ 

carefully  kept  in  mind.  The  function  .  r  represents  tne 
torque  exerted  by  sun  and  moon.  It  is  the  same  for  a  rigid 
and  a  nonrigid  earth.  The  function  . is  caused  by  elast 
deformation;  it  therefore  depends  on  the  Love  number  k  and 
is  zero  for  a  rigid  body,  in  accordance  with  ; 4  -  1 2  )  .  Thus 
„  represents  the  "direct  tidal  efte.t"  in;  ....  rhe  "  in¬ 
direct  tidal  effect",  and  the  sun; 


expresses  the  total  effect  of  sun  and  moon. 

Solution  of  the  basic  equation.  in  view  of  (5-2  1  a n 0 
(5-21),  the  basic  equation  (5-1)  may  be  writ1-. en 


rfi  -  i  7  E  (  m 


RD  ' 


Using  (4-3)  and  (4-5),  this  can  be  brought  into  the  form 


rfi  -  i  t  m 
c 


-  1  - 


1  + 


k  -E  ’  T 


(5-23 


,  k  !• 

i+r  - 


( 1  + 


IB  e' 
1 


The  solution  of  this  equation  is  given  by 

,  k  ^ 

it)  =  m  ~  e  1  c  ‘  +  i  -  — — —  — —  B  e  ~ 

■  -j  +  i  ,  k  i 

T  l  C  l+i-  — ^  J 

as  can  be  verified  by  substitution.  It  contains  two  constants, 
which  are  the  real  and  the  imaginary  part  of  the  complex  con¬ 
stant  in  j  ,  and  is  thus  the  general  solution  of  (5-23). 

The  solution  (5-24)  consists  of  the  solution  of  tne 
homogeneous  equation  (4-4),  co r res po nd i n g  to  free  motion 
without  external  forces,  and  a  term  representing  the  effect 
of  lunisolar  perturbations. 
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POLAR  MOTION  :  OTHER  AXES 


6  . 


Eq.  (5-24)  describes  the  movement  of  the  instanta- 
neous  rotation  axis,  which  is  c h a rac te r i 2 e d  by  the  complex 
number  m  =  m ;  +  im2  where  in j  and  in  are  the  equatorial 
components  of  the  unit  vector  of  the  rotation  axis.  Similarly, 
the  angular  momentum  axis  is  characterized  by  the  complex 
number  h  given  by  (3-25),  and  the  figure  axis  is  described 
by  the  complex  number  f  given  by  (3-27). 

Eq.  (3-27)  gives  the  figure  axis: 


where  c  =  C; ;  +  iCj  j  is  a  complex  combination  of  proaucts 
of  inertia.  By  (3-22),  (3-23),  and  (3-25)  we  have 

h=^m+^c.  (6-2) 

The  combination  of  (6-1)  and  (6-2)  finally  gives 
,  A  C-A  c 

h  =  m  +  — £—  f  (  6  -  3  ) 

for  the  angular  momentum  axis. 

For  the  resulting  formulas  we  shall  recall  the  notations 
for  a  few  constants  and  shall  introduce  a  new  one  :  the 
secular  Love  number  (4-2) 
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the  Eulerian  frequency  (3-13) 


C-A 

~r  - 


( 6-5  ) 


the  Chandlerian  frequency  (4-5) 


1-* 


(6-6; 


the  nutational  frequency  (5-20) 


3  3 

(the  name  will  become  clear  in  sec.  7),  and  the  factor 


(6-7) 


i  +  i 

K 


i  c  k  _ 


!  1  +  f 

i  k 


T< 


(  6  -  8  ) 


Using  these  notations,  we  can  represent  the  polar  motion 
p  of  various  axes.  From  (5-24)  we  get 


i  v-  .C-A  ,■  D  - 1.  ( -jj 
m  =  rn  e  +  i  — '-.Be  i 

A  -.11 


6-9  1 


for  the  polar  motion  of  the  rotation  axis  R  .  The  f i gu re 
axis  F  is  obtained  from  (6-1),  where  c  is  the  sum  of 
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(4-1)  for  rotational  deformation  and  (5-17)  for  tidal  defor¬ 
mation  : 


P„  = 


f 


k 

k 


m,  e 


1 


c  L 


(6-10) 


Finally,  these  two  equations  are  linearly  combined  by  (5-3)  to 
getthe  polar  motion  of  the  angular  momentum  axis: 


PH  =  h 


.  i  a  -  t 

h  ne  c  +i 


C-A 

C 


1  -  — 
*•  L, 


—  B  e 


c  +  -  ) 


7  1  +  V - 


(6-11) 


with 


n 


k  1  m 
kT  !  0 

S 


(6-12} 


These  formulas  are  illustrated  by  Fig.  6.1,  which  shows 
the  plane  tangent  to  the  terrestrial  sphere  at  the  point  0 
representing  the  z-axis  and,  at  the  same  time,  a  mean  position 
of  all  three  axes.  The  other  points  designate  the  rotation  axis 
R  ,  the  angular  momentum  axis  H  ,  the  figure  axis  F  ,  and 
their  force-free  counterparts  R  , ,  H0  and  F  ,,  as  given  by 
the  fi-st  term  on  the  right-hand  side  of  (6-9),  (6-11),  and 
(6-10'  .  In  order  to  get  a  feeling  for  the  orders  of  magnitude, 
we  note  that 
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s 


C-A  .  C-A  . 

~T  =  ~J~  = 


0.003 


(5-13; 


Let  us  first  consider  the  free  motion  only.  The  points 
R 3  and  H 0  describe  concentric  circles  whose  radius  is  on 
the  order  of  0.2“  which  corresponds  to  6m  .  Both  points 
are  very  close  to  each  other  since  h0  =  by  (6-12)  and 

(6-13);  there  is  HoR0  =  2cm.  The  point  0  corresponds  to 
the  figure  axis  in  the  undisturbed  case.  Since  the  axis  of 
rotation  does  not  coincide  with  the  figure  axis,  the  rotation 
of  the  earth  produces  a  nonsymmetnc  deformation  (sec.  4)  which 
causes  the  axis  of  maximum  inertia  to  shift  to  F;  .  The  first 
term  on  the  right-hand  side  of  (6-10)  shows  that  OF  4  2 in 
since  k  /  k  s  =  0.3.  The  prints  Rq  ,  H0,  and  Fi  lie  on  the 
same  radius  and  slowly  rotate  together  around  0;  the  perioc 
is  the  Chandler  period  of  about  430  days. 

So  much  for  the  free  motion.  The  attraction  of  sun  and 
moon  causes  forced  motions  which  are  represented  Dy  the  second 
term  on  the  right-hand  sides  of  (6-9),  (6-10),  and  (6-11).  The 
instantaneous  pole  of  rotation,  R  ,  describes  a  near-circular 
closed  curve  around  Rc ;  we  have  R.R  =  60  cm.  A  similar 
curve  is  described  by  the  angular  momentum  pole  H  around  H ■ ; 
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since  1  -  k  /'  k  =  0.7,  we  have  H^H  =  40  cm.  Especially 
remarkable  is  the  motion  of  the  pole  F  of  the  figure  axis: 
it  describes  a  quasi-circular-  motion  around  Fo  whose  radius 
is  by  (k/k  )A/(C-A)  =  100  times  larger  than  ORq  ,  namely 

about  60  in  !  From  (6-3)  it  is  clear  that  R,  H,  and  F  lie 

on  a  straight  line. 

Since  ^  =  .. ,  the  period  of  these  forced  motions  is 
on  tne  order  of  1  day;  we  therefore  speak  of  diurnal  polar 
motions . 

We  also  point  out  the  evident  fact  that  the  motion  of 
the  various  axis  for  the  case  of  a  rigid  earth  can  be  obtained 
by  putting  the  Love  number  k  =  0  in  these  formulas.  The  main 
change  in  Fig.  6.1  is  that,  for  a  rigid  earth,  the  points  F 

and  F-  will  coincide  with  the  origin  0  and  the  radii  R:R 

and  H '  H  will  be  almost  equal. 

A  final  remark  is  in  order.  The  actual  free  polar  motion 
is  much  more  complicated  than  the  simple  circular  model  conside 
ed  here,  for  a  variety  of  reasons,  not  all  of  which  are  well 
understood.  The  free  motion  cannot,  therefore,be  adequately  pre 
dieted  by  an  analytical  model  and  can  only  be  determined  by 
observation  (International  Latitude  Service,  International  Pola 
Motion  Service,  Doppler,  Laser,VLBI).  On  the  other  hand  it 
appears  that  the  lunisolar  (forced)  motion  of  the  pole  can  oe 
pred i c  ted  well. 
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Tni:  equation  is  differentiated  with  respect  to  time,  we  trea 
H  as  a  constant  since,  to  a  high  accuracy, 

H  -C  =  const.,  (  "  -  3 


c  f .  tne  third  equation  of  (2-11).  The  first 
of  (7-1)  are  then 
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’hen,  two  symmetric  frequencies  -  .  and  give  tne  same 

nutations!  frequency  ^  (only  tne  sign  of  the  coefficients 

J 

will  be  different). 

It  is  clear  that  the  sidereal  frequency  ..  =  produce 

precession.  The  principal  nutation  component  comes  froi!i  an 
that  corresponds  to  the  motion  o*  the  ’unar  node.  It  has  a 
period  of  13.65  yean  and  amplitudes  in  longitude  (  i  .  e  .  ,  ir.  • 
o f  17.2"  a  ri  d  in  obliquity  (  i  .  e  .  ,  in  )  of  9.2"  .  Tne  1 1  d  a  ’ 
effect  of  this  ,u.t;on  is  insignificant.  Put  in  view  of  the  smal 
ness  of  j. ...  ,  the  corresponding  nutations!  coefficient,  proper 

i  o  n  a  1  to  3  _  /  .j  b  y  '7-14;.  is  grea  1 1  v  magnified.  0  the  r  p  e  >■  i  o  d 


in  precession 


nutation  figuring  in  (7-14)  are  simply  related  to  tr.e  coeffi¬ 
cients  of  tne  tidal  potential,  namely  through  (5-14;.  The 
fact  that  precession  (as  well  as  nutation;  is  proportional 
to  the  cons  tant 


makes  it  possible  to  determine  this  constant.  It  is  called 
dynamical  ellipticity  and  is  of  basic  importance  for  physical 
geodesy  ( c  f .  Heiskanen  and  Moritz,  1967,  p.339). 

The  relation  between  precession,  nutation  and  tidal 
potential  has  been  studied  in  particular  detail  by  Melchior 
(1971,  1978). 

Finally  we  point  out  that  precession  and  nutation  of 
the  angular  momentum  axis  depend  only  on  the  lunisolar  torque 
and  not  on  the  Love  number  k  .  Therefore,  the  formulas  (7-14' 
are  the  same  for  a  rigid  and  an  elastic  earth;  cf.  also 
(Fedorov,  1963,  p.  16). 


3.  PRECESSION  AND  NUTATION  :  OTHER  AXES 


1 

I 


Of  the  three  axes  considered  in  sec. 6  :  the  instanta¬ 
neous  rotation  axis,  the  angular  momentum  axis,  and  the  figure 
axis,  we  have  treated  the  spatial  motion  of  the  angular  momen¬ 
tum  axis  in  the  preceding  section.  The  spatial  motion  of  the 
remaining  two  axes  will  be  studied  now.  As  we  have  seen  in 
sec. 6  ,  however,  the  instantaneous  figure  axis  performs,  in 
the  case  of  an  elastic  earth,  such  a  large  daily  motion  with 
respect  to  the  earth's  body  (about  60  m )  that  it  is  of  little 
practical  use.  Much  more  useful  is  the  z-axis  which  corresponds 
to  the  figure  axis  of  an  undeformed  earth  ana,  for  a  rigid  eart 
coincides  with  the  figure  axis.  So  we  shall  study  the  motion 
of  the  z-axis  and  of  the  i n s ta n ta neo us  rotation  axis,  as  well 
as  the  so-called  celestial  pole. 

It  will  turn  out  that  the  precession  of  different  axes 
is  the  same,  and  the  nutation  nearly  so.  The  very  small  differ¬ 
ences  between  the  nutation  of  the  angular  momentum  axis  and 
that  of  other  axes  are  sometimes  called  Oppolzer  terms  (Woolard 
1953;  Kinoshita,  1977). 

Motion  of  the  z-axis.  We  again  employ  our  usual  two 
coordinate  systems  :  a  "space-fixed"  inertial  coordinate  system 
X  ;  X  :  X  3  as  specified  in  the  beginning  of  sec. 7  (cf.  Fig.  7.1 
and  an  "earth-fixed"  system  x;x:x;  introduced  in  sec. 3.  In 
the  latter  system,  the  axis  x  ,  =  z  corresponds  to  tne  figure- 
axis  of  the  undeformed  earth,  and  it  is  a  L i ouv i 1 1 e  axis,  the 


equation  of  motion  taking  the  simple  form  (3-8).  The  z-axis 
is,  tnerefore,  also  called  a  mean  Tisserand  figure  axis.  The 
X  •_  X  plane  is  the  ecliptic,  and  the  X  <  axis  represents  the 
vernal  equinox  (both  at  a  fixed  epoch).  The  x : x ?  plane  re¬ 
presents  the  equator  (more  precisely,  the  mean  equator  of  fi¬ 
gure),  and  the  Xj  axis  corresponds  to  the  Greenwich  meridian 
(more  precisely,  to  a  conventionally  assumed  fixed  direction 
close  to  the  Greenwich  meridian). 

Since  we  are  concerned  only  with  directions  and  rotations 
the  origins  are  of  no  interest  here;  we  can  for  the  present  pur¬ 
pose  consider  both  systems  geocentric. 

The  relative  orientation  of  the  x-xx.  system  relative 
to  the  X,X .Xv  system  can  be  given  by  the  three  Euler  angles 


FIGURE  8.1.  T  h  e  basic  Euler  a  n  q 1 e  s  . 
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;  ,  9  ,  v  defined  as  in  Fig.  8.1.  The  quantities  .  ,  ■  ,  are 
similar,  but  not  identical,  to  the  angles  ,  :  ,  of  Fig. 

ti  ;i  i'. 

7.1  :  now  the  pole  is  the  z-axis  and  not  the  angular  momen¬ 
tum  axis.  We  get  the  system  x ; x  x  by  rotating  the  system 
X;X;X'3  first  about  the  3-axis  by  the  angle  .  ,  until  it 
coincides  with  the  node  N  (which  is,  of  course,  not  exactly 
the  same  point  as  N  in  Fig.  7.1),  then  about  the  nooal  axis 
by  the  angle  -•  ,  and  finally  aoout  the  3  -  a  x  l  s  by  the  angle 

The  angle  u  is  the  longitude  of  tne  node,  repre¬ 

sents  the  obliquity  of  the  ecliptic,  and  :  is  an  angle 
that  measures  the  rotation  of  the  earth.  These  Euler  ancles 
are  frequently  used  in  physics  and  astronomy;  note,  however, 
that  almost  every  author  uses  a  different  definition  of  them. 
We  follow  (Plummer,  1918)  and  (McClure,  1973). 

The  components  of  the  rotation  vector  ..  (wnic.n  does 
not  in  general  coincide  with  the  x  .  =  z  axis!)  are  oenotec 
by  .  Compared  to  (3-11),  (3-12)  and  (3-13)  we  nave 

w  i  =  m  i  ,  ^2  =  ,  v  j  =  .  {  1  +  m  .  )  .  ■  8  -  i  ) 

They  are  connected  to  the  time  derivatives  of  the  Euler  angle 
by  Euler's  well-known  kinematical  equations: 

j  =  -  9  c  o  s  o  -  q  s  l  n  c  s  i  n  .*>  , 

cj  2  =  9  s  i  n  i  -  o  s  i  n  9  c  o  s  ?  ,  (3-2  ; 

=  vCO  S 9  +  p  , 


"  3 


c  f .  (Synge,  1960,  p . 2 8 )  ,  but  note  the  different  definition 
of  Euler  angles. 

The  complex  combination  of  the  first  two  equations 

gives 


« .  +  =  -e  1  *' (  9  +  i/Sinr  )  ,  (  3  -  3  ) 

w  n  e  n  c  e 

-  +  i * s  i  n 1  =  - e  1  f  „ ;  +  i  ) 

=  -  e  1  f ..  (  m  •.  +  i  m  _  )  , 

by  (3-1)  ,  or  briefly,  using  (3-15) 

•  +  i  1  s  i  n  •  =  -..me  1  *  .  (  8-4  ) 

This  equation  is  of  basic  importance:  it  relates  precession 
and  nutation,  as  expressed  by  the  Euler  angles  .  and  -  ,  t 
the  polar  motion  of  the  rotation  axis,  m  . 

The  angle  ?  measures  the  earth's  rotation.  Since  tne 
changes  of  ••  and  ;  are  very  small,  it  is  a  sufficient  ap¬ 
proximation  to  put,  similarly  to  (7-10), 

=  i-t  .  (3-5) 

Then  we  can  substitute  m  from  (6-9)  into  (8-4)  to 


obtain 


+  i  v  sin- 


..m  i  e  1 


.  B  e  ’  1  '  '  ‘  i 


This  equation  is  very  similar  to  (7-13)  ana  is  integrate!  in 
the  same  way,  with  the  result 


+  i  a v s  i  n -  -  i  — +r—  in  - e 


B  t 


C-A  -  „ 

IT  '  IT " 


where,  in  analogy  to  (7-14),  ,\  and  t .  are  the  d  i  f  f  e  >-e  n  o  e 
between  the  instantaneous  values  of  ■■  jnd  ,  and  some  re¬ 
ference  values. 

This  equation  is  very  similar  to  (7-14).  in  fact,  the 
precession  term  is  the  same,  and  the  nutation  terms  only  dif 
fer  by  the  factor  <  which  is  very  close  to  unity  by  >  6-3 
and  (6-13).  An  important  difference  is  the  presence  of  tne  f 
term  on  the  right-hand  side  of  (8-7),  which  originates  from 
free  polar  motion. 

This  free  motion  term  has  a  nearly  diurnal  frequency, 
since  .3  -  :  is  very  close  to  the  sidereal  frequency.  The 

forced  motion  (lunisolar)  term,  on  the  other  hand,  is  long- 
periodic,  the  principal  periods  being  18.66  years,  a  year, 


between  the  Euler  angles  for  the  rotation  axis  and  those  for 
the  z-axis  as  considered  in  (8-7). 

Fig.  8.3  shows  that  v  and  (  -  '  vs  i  n  •; )  are  related 
to  the  polar  motion  components  m,  and  m:  by  a  plane  ro¬ 
tation,  which  is  best  written  in  complex  form: 

•  +  i  j  s  i  n  =  -  i  me  1  ’  =  -ime1"t  (8-3 

where  m  =  m .  +  im2  as  usual.  Again  we  substitute  (6-9)  and 
obtain 


i  ( -  i  +  i  >  t 

+  i  ■  ,.  s  i  n  ••  =  -  i  m  e  l-  + 


C  -  A  -  p  -  i  (  Z  uj  t  li  . 
TT  ■■  ■  jBje  J  ’ 


(S-l 


This  is  added  to  (3-7)  and  gives 
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The  equation  is  very  similar  to  '  >’-14  j  and  (1-7;,  A  q  a  ’  n  , 
the  precession  term  is  the  same,  and  the  nutation  c oe r f i c ; en t j 
nearly  so.  Again,  there  is  a  term  due  to  free  polar  motion  toe 
first  on  the  right- h and  side)  which,  however,  is  much  smaller 
than  the  free  motion  term  in  (3-7). 

The  celestial  pole.  Tree  polar  motion  cannot  be  ad¬ 
equately  modeled  in  a  simple  mathematical  way;  the  circle 

m o e  * ' 

around  0  with  radius  OR'  (Fig.  6.1)  is  valid  only  for  an 
ideally  elastic  earth.  In  reality,  the  polar  motion  curve  is 
rather  irreqular  and  can  only  be  determined  empirically;  cf. 
(Mueller,  1969,  p.33). 

On  the  other  hand,  lunisolar  effects  can  be  predicted 
well,  both  in  polar  motion  and  in  nutation.  It  is,  therefore, 
appropriate  to  refer  calculation  of  precession  and  nutation  to 
an  axis  which  is  not  affected  by  free  polar  motion.  This  ex¬ 
cludes  the  rotation  axis  (3-11)  and  the  z-axis  (3-7),  but  lea¬ 
ves  the  angular  momentum  axis  which  does  not  contain  a  free 
motion  term  in  nutation  by  (7-14). 

An  inspection  of  Fig.  6.1.  shows,  however,  that  the  in¬ 
stantaneous  angular  momentum  axis  H  shows  a  forced  nearly 
diurnal  polar  motion  of  radius  40  cm  around  H  >  .  If  we  wisn 

to  have  an  axis  which  not  only  is  a  suitable  reference  for  pre¬ 
cession  and  nutation,  but  also  is  relatively  stable  with  respect 


mi*[i  rAi.t 


- -■-*»-** 


Thus,  by  '  S  -  1 2  )  ,  the  po  !  a  r  mo  T  i  jn  of  ,1  cun  *  a  i  " 

*  o  reed  part  and,  by  '8-  1  *2  :  ,  it;  n  u  t  .i  t  i  >>  n  ■,m;  t  a  i  r.  •,  >,<,  f  *- *r  <- 

oar”:..  3ot.n  missinu  carti  would  hav>  had  a  near!  /  J  ;  j  r  r; 1 
period  as  we  have  pointed  o  u  t  in  above.  mo  t  i  on  ,  a  >• 

thus  optimally  smooth. 

Furthermore,  C  =  H  shares  with  the  ang  u  1  a  r  ;n  c  men 

turn  axis  H  the  property  that  its  nutation  i  a  n  a  i  y  t  r  •.  a  1 

predictable  since  it  does  not  contain  a  tree  te»*m  w ren 
would  be  accessible  only  afterwards  by  jbse**va  1 1  •  >  n  . 

c or  these  and  other'  reasons  ■:  L e  i  c  x  .  1878;  L e  i  c *  arc 

Mueller,  1979),  it  appears  that  C  is  try  best  candidate 
an  appropriate  definition  of  a  celestial  reference  pu'e. 
its  General  Assembly  in  Montreal  in  le  7  a.  the  Intemati  ; 
Astronomical  Union  has,  in  fact,  adopted  U  as  toe  u  •  *  i 
celestial  n o 1 e  of  reference. 

We  finally  note  that,  for'  m  -  0  ,  t n e  expressions 
'8-7  and  (8-15)  coincide.  Eq .  (6-10)  shows  that  for  a  >’ i 
body  (  k  =  0  ,  also  f  --  0  in  this  case,  so  t  n  a  t  t  n  e  f  i  gu  re- 
ax:  s  remains  at  the  origin  of  polan  motion.  Thus,  for  the 
of  a  rinid  earth  and  in  the  absence  of  noiar  mot' on  out 
generally!),  the  celestial  Dole  coincides  wifi  the  o ;  - 
o  f  the  f  i  ou  re  axis.  There  fore  ,  the  d  i  ret  t  i  on  o  r  Me  e  1  •  ■  . 
pole  nas  also  been  called  "figure  a  x  1  .  T,r  unfortunate 

terminology  (Atkinson,  1976)  has  c  a  u  s e d  c  o n  -  i J e r a  P i  *  •  c  .  n  * 


as  des c  r i 

bed 

by 

Mue 

1  1  e  r  (  1 9 80) . 

I  f 

the 

ea  r  th 

is  elastic  (  e 

v  e  n  wit >' 

a  1  i  u  i  d 

dc  1  a  r  tio  t 

i  on 

i  s 

abs 

e n  t ,  then  we 

c.  an  o  n  1  / 

,  a  v  that 

w ; t  h  the 

o  r  i 

a  i  n 

0 

in  Fin.  6.1, 

w  h  i  .;  h  i  1 

.  seen  to 

nutation  has  some  affinity  to  the  present  m-' nod  ,  esuec  j  .t  I  1 
in  tne  distinction  between  annular  monion  t  urn  a  /.  i  •  ,  figure- 
axis,  and  rotation  axis;  this  relation  ha  been  pointed  .j  s  t 
by  McClure  (19  7  3,  sec.8).  It  is  restricted  to  a  rigid  e-jrtri 
and  presents  a  less  systematic  and  complete  treatment  o T  a ’ 
aspects,  but  numerical  details  suen  as  the  neglect  of  terms 
are  more  carefully  considered  (although  several  authors  h  a  v 
pointed  out  minor  en  ors  in  W  o  o  1  a  r  d  '  s  treatment),  w  o  o  1  a  r  d  1  s 
theory  has  b«en  a  breakthrough  and  a  classic.  It  has  served 
for  decades  as  the  official  reference  for  precession  and 
nutation  of  the  International  Astronomical  Union.  Now, 
however,  it  has  been  superseded  by  Kinosnita's  theory  for 
a  >'igid  earth  and  by  other  theories  for  elastic  and  liquid 
co  re  models. 


k:noshita's  theory  of  precession  and  nutation. 


The  most  elegant  formulation  of  analytical  mechanics 
is  Hamilton's  theory.  It.  can  be  summarized  as  follows,  cf. 
(Arnold,  19  73)  or  (Synge,  1960).  Let  a  conservative  mecnar.i- 
c a  1  system  be  descrioed  by  n  independent  variables  q. 
which  are  called  "generalized  coordinates",  let  its  kinetic 
energy  be  T  and  its  potential  energy  U  .  Define  the  "ge¬ 
neralized  impulses"  p  by 


w  n  o  r  a  g  -  dp  /  d  t  ,  and  define  the  Hamiltonian  function 
b  t 

9-2 ; 

Then  the  equations  of  motion  for 
consideration  are  Hamilton's  canon 


•it  "  •;> 


eo ual  to  the  total  energy, 
the  ivnamical  system  under 
. a  I  equations 


dp.. 

■It  ■<!. 


(9-3) 


* 

i 

The  quantities  p  and  q,  are  called  canonical  v r  :  a  h  1  ,  . 

Because  of  its  formal  simplicity,  Hamilton's  ineory 
possesses  considerable  theoretical  advantaqes.  canonical 
variables  are,  therefore,  frequently  used  in  celestial  me¬ 
chanics;  cf.  (Brouwer  and  C  lenience,  1961,  chapter  XVII,. 

They  have  been  introduced  into  the  theory  of  precession  and 
nutation  by  H.  Andoyer  in  1911;  cf.  (Andoyer,  1  9  23,  192^]. 

Recently,  Kinoshita  (1977)  has  used  Andoyer  variables  to  de¬ 
rive  the  most  accurate  theory  of  precession  and  nutation 
available  for  a  rigid  earth.  In  this  section  we  shall  present 
the  theoretical  foundation  of  Kinoshita's  theory  which  is 
very  simple  whereas  the  details  are  enormously  comp  1 i ca ted . 

Andover's  variables  are  denoted  by 

Q  1  =  1  P ;  -  L  , 

q :  =  g  ,  P =  G  ,  (9-41 

q  5  =  h  ,  P  ;  =  H  , 

so  that  coordinates  and  correspondinn  impulses  are  denoted  'ey 
the  same  letter.  Their  definition  is  as  explained  by  means 
of  Fig.  9.1. 

We  use  the  two  systems  XX.X,  (inertial)  and  x:x  x; 
(earth-fixed)  as  before  (cf.  sec. 8).  The  ecliptic  (fixed  at 
an  epoch  t  ^  )  corresponds  to  the  X  L  X  .  plane,  and  the  (instan¬ 
taneous)  equator  to  the  x  ;  x  ■  plane.  The  point  x.  denotes 
the  figure  axis  which,  in  fact,  for  a  rigid  body  coincides  with 
the  z  =  axis.  The  point  M  denotes  the  angular  momentum 
vector  M  .  Note  the  change  in  notation  for  this  vector  which 


FIGURE  9.1.  The  A  n  d  o  y  e  r  variables  n  ,  g,  1  .  .he 

point  M  denotes  the  angular  momentum 
vec  tor . 


elsewhere  in  this  report  is  denoted  by  H  .  This  cnange  is 
necessary  since  we  wish  to  retain  Kinoshita's  notation  (  9  -  A 
where  H  is  a  canonical  v  a  r i a  b 1 e  .  The  equation  of  M  cor 
resounds  to  a  plane  which  is  normal  to  the  angular  momentum 

vec  t  o  r  M  . 

Fiq.  9.1  can  be  considered  a  superposition  of  Figs.  7 
and  .  1  ,  with  the  important  difference  that  now  we  look,  s 
tr,  s  a  r  .  at  the  back  of  the  unit  sphere  :  the  nodes  N  an 
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Q  shown  in  Fig.  9.1  differ  from  the  corresponding  nodes  in 
the  previous  figures  by  180 

Apart  from  this  difference,  the  Andoyer  variable  h 
corresponds  to  ;■  in  Fia.  7.1  and  the  new  variable  1  is 

seen  essentially  to  be  the  argument  of  polar  motion  (  ■  t  *  ■ 
in  eq.  (2-16)).  The  sum  g  +  1  is  nearly  equal  to  the  Euler 
angle  :  +  180°.  The  angles  I  and  I  in  Fig.  9.1  are 
equivalent  to  e  and  e  in  the  previous  figures. 

Kinoshita  also  introduces  the  spherical  distances 

X,N  =  hf  ,  N  x i  =  p  ,  (9-5) 

correspondi  ng  ,  in  Fig.  8.1.  ,  to  -p  and  :■  ,  re  s  pec  t  i  ve  1  y  . 

Hence  h  ,  I e ,  $  are  the  Euler  angles  (apart  from  180  )  in 
Kinoshita's  notation. 

Having  thus  introduced  the  canonical  coordinates  1,  g 
h,  we  can  easily  find  the  correspondi ng  canonical  impulses 
L,  G,  H  by  (9-1).  First  we  have  to  find  the  kinetic  energy 
T  .  From  classical  mechanics  (Arnold,  1978,  p.  137)  we  know 
the  relation 


where  r  is  the  inertia  tensor  (2-10),  j  is  the  rotation 

T 

vector,  and  ^ '  its  transpose.  The  vector  ^  can  now  be  writ 
as  follows: 


j  =  hel*ge  <-  1  e ,  . 

— h  —  q  —  L 


ten 


(  9  -  7  ) 


6  b 


This  means  that  the  total  rotation  is  split  up  into  a  rotation 
about  t  n  e  X .  axis  with  speed  n  =  d  h  /  d  t  ,  into  a  rotation 
about  the  angular  momentum  axis  witn  speed  g  ,  and  a  rotation 
about  the  X.  axis  with  speed  1  .  The  vectors  e.  ,  e  „  , 
and  e,  are  the  unit  vectors  of  the  directions  around  whicn 
the  rotations  are  performed:  e  is  the  unit  vector  of  the 

—  n 

X-.  axis  (the  angle  h  is  counted  in  the  plane  normal  to  tnis 
axis),  is  the  unit  vector  of  the  angular  momentum  vector 

M  (the  angle  g  is  counted  in  the  plane  normal  to  M  ),  and 
e,  is  the  unit  vector  of  the  x.  axis  (for  a  similar  reason) 
The  differentiation  of  (9-6)  with  respect  to  1  yields 


=  M  •  e,  ,  (9-8) 

because  of  the  symmetry  of  C  ,  using  (2-5)  with  M  instead 
of  H  and  (9-7),  and  denotinq  the  inner  product  by  a  dot  as 
usual.  Other  d  e  r  ■  v  a  t  i  v  e  s  are  obtained  in  the  same  way.  The 
canonical  impulses  are  now  given  by  (9-1)  with  (9-4)  : 

,  T 

L  -  -  =  M  •  e,  , 

*  i 

o  -  -v  --  M  •  ea  ,  (9-91 

-  n 

-•  -  M  •  e 

•  h 


n 


'bus  ,  L  is  the 


component  of  M  ,  G  is  the  magnitu- 


of  the  annular  momentum  vector  itself  since  :n..- 


unit  vector 


and  M  have  the  same  direction 


is  the  X 3  component  of  M  . 

Briefly  we  may  say  that  a  canonical  impulse  'L,  G, 
or  H)  is  the  component  of  M  normal  to  the  plane  along  wnicn 
the  corresponding  canonical  coordinate  ( 1 ,  g,  or  hi  is  counted 
Using  the  angles  I  and  J  in  Fig.  9.1  ,  we  may 

w  r  i  t  e 


G  =  1  M1  , 


L  =  GcosJ  , 


H  =  GcosI  .  (9-10; 


If  we  know  the  canonical  variables  (1,  g ,  h,  L,  G ,  H), 
then  the  Euler  angles  can  easily  be  computed:  by  (9-10'  we 
have 


cos  I  =  £  , 


cos  J  = 


(9-11 


the  solution  of  the  spherical  triangle  NPQ  then  gives 
QN  and  NP  ,  and  finally 

h,  =  X  i  N  =  h  +  QN 
j=Nx-  =  N  P  +  1  . 


We  now  need  an  expression  of  the  kinetic  energy  T  in 
terms  of  canonical  variables.  In  the  body -fixed  system  x  ;  x ;  x 
formed  by  the  principal  axis  of  inertia,  the  expression  (9-6) 
becomes,  using  (2-10)  and  (2-11)  , 


:  -  *  vi  •  .  t  ■  M  -  *'■  v 

•  .  *•’  .  i  .  ■  q  .  •  , 

2  '  A  S  C  '  ’  '  -- 

we  do  n  o  t  l>  r  e  $  u  3  p  o  so  r-otat  ionjl  symmetry  so  tna  t  toe  p  r  i  r.  c  i  - 
pal  e  g  a  a  to  r  i  a  i  mo  me  n  r  s  of  inertia,  A  and  E  ,  may  be  di  f ‘er¬ 
ect  . 

Sow,  M  .  -•  L  by  (9-9  )  .  Further,  the  projection  of 
o n to  the  x • /  o'  a n e  i  s 

.  M;  +  M  .  G  -  L  3.  1  _ 

him  ‘nr  ms  with  the  x  _  axis  the  annle  270  -  1  ,  so 

t  n  a  t 


M  -  -  G  ■  -  L  ■  s  i  n  1  , 
M  -  -  .  G  ~  -  L •  cos  i 
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^ e n ce ,  9-13  )  beco me s 


sin-  1 


cos  •  1 


(G-  -Li 


2  C 


(9-16) 


which  expresses  the  kinetic  energy  in  terms  of  the  canonical 
v  a  "*  i  a  h  1  n  s  . 


Now  we  may  write  the  Hamiltonian  equation’',  '  )  -  3  i  r  o  r 
our  problem: 


dl 

s - 

dJL 

3T  " 

'•T  ’ 

d  t 

•c~T 

•g  ... 

~  y 

dG 

i  ’• 

Jt  " 
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d  t 

•  a 

9h 

"  Tv 

dH 

n 

"7h  ’ 

dT  r  ' 

sTT 

3y  (9-2)  we  have 
.V  =  T  e  U  , 

where  I  is  given  by  (9-16)  and  the  potential  energy  U  comes 
from  the  attraction  of  sun  and  moon. 

Equations  (9-17)  may  be  solved  by  standard  perturbation 
methods  of  celestial  mechanics  (Brouwer  and  Clemence,  1961, 
chapter  XVII),  usinq,  as  a  first  a pp ro x  i  ma t  i  on ,  free  motion  with 
U  =  0  .  This  free  motion  is  discussed  in  detail,  e.g.  in 
(Arkhanqelsky,  1977).  Andoyer  (1926)  used  simple  variation  of 
constants,  whereas  Kinoshita  employs  a  considerably  more  sophi¬ 
sticated  technique  (method  of  Hori,  cf.  (Schneider,  1979,  p.AQ9i 
We  cannot  go  into  the  details  of  this  solution  of  the 
forced  motion,  which  are  enormously  envolved,  and  refer  the 
reader  to  (Kinoshita,  1977).  The  results  are  series  for'  precessi 
and  nutation  of  the  anqular  momentum  axis,  the  figure  axis 
(for  a  riqid  body  coinciding  with  the  z  =  x:  axis),  and  the 


rotation  axis  which  have  the  form  of  tne  series  derived  in 
secs.  7  ana  3  .  In  fact,  ;.h  =  ,\v..  and  '.I  =  are 

h  r. 

nutation  in  lonqiti.de  and  obliquity  for  tne  angular  momentum 
axis,  and  similarity  for  the  otner  axes. 

Kinoshita's  results  are  accurate  to  0.0001"  ,  corre¬ 
spond!  nc  to  3  mm  in  position.  They  represent  the  most 
complete  and  precise  theory  of  precession  and  nutation  a  v a i 1  a : 
for  a  rigid  earth,  especially  because  he  developes  a  very 
accurate  expression  for  the  lunisolar  potential,  and  are  tnus 
a  progress  with  respect  to  'VJoolard,  1953).  Kinosnita's  t  h  e  o  r  > 
is  also  somewhat  more  accurate  than  the  method  described  in 
secs.  7  and  8,  but  this  latter  met  nod  is  valid  also  for  an 
elastic  earth  whereas  Kinoshita's  method  is  restricted  to  a 
rigid  body. 

The  best  way  to  compute  nutations  for  an  elastic  earth 
seems  to  apply  the  formulas  of  sec.  8  to  compute  differences 

"  alasii;  r i  i  i  i 


by  applying  these  formulas  first  for  the  actual  Love  number 
k  and  then  for  k  =  0  (which  gives  rigid  body  results). 
Since  these  differences  are  very  small,  they  can  be  computed 
quite  precisely  by  the  theory  of  sec.  8  .  These  differences 
are  then  added  to  Kinoshita's  rigid-body  results  to  dive 


A  similar  procedure  is  advocated  by  Wahr  (1979,19d0 
to  take  into  account  effects  of  the  liquid  cor-e  .  In  fact, 
such  effects  are  considerably  larger  and  more  important  tn 
effects  of  elasticity.  Therefore,  the  remaining  part  of  tn 
report  will  be  devoted  to  the  influence  of  the  liquid  core 
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'hp  p  I’isiderar.  ion  of  the  free  rotation  of  a  rigid  body 
i :  .3  n  ••  ’  oen  va  1  •,<-  problem  is  of  basic  theoretical  importance 
m  •  '  1  >erve  as  a  preparation  for  a  detailed  understanding 
or  ’Maid  cor-  effects  to  be  treated  later  in  this  report. 

We  snail  restrict  ourselves  to  a  rotational ly  symmetric  earth 
'  .*  r  i  :  n  tne  principal  equatorial  moments  of  inertia  are  eq  ua 


A  >■  3  ;  (10-1 

^  general  case  A  =  B  can  be  treated  in  a  similar  way. 

eigenvalues  for  Euler’s  Equations.  By  (2-13)  we  have 


A  1  ;  f  f  C  -  A  )  ■.,  „  :  0  , 

A  : .  -  ; c -  A ) . , ,  =  0  ,  (10-2 

C  :  .  =  0  . 

Here  C  is  the  maximum  (polar)  moment  of  inertia,  and 
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where  represents  a  constant  value  for  the  average  s u e e d 

rotation  and  ~ ^ :  are  very  small  as  compared  to 

that  their  s  q  u  a  r e  s  ana  higher- order  powers  can  be  neglected 
Thus  we  may  put  _  .  3  i n  the  first  two  equations 
(10-2),  so  that  Euler’s  equations  become  linear: 


=  0  , 

=  0  , 

=  0  , 


;  io 


where 


C  -  A 


is  the  Euler  frequency  (2-14)  .  'we  are  just  looking  at  the 
problem  of  sec.  2  from  a  somewhat  different  angle.) 

As  before,  it  will  be  convenient  to  use  complex  ouan 
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T  h  i  s  is  a  very  simple  quadratic  equation  for  ,  o  r'  w  h :  c  r 
the  solutions  obvious1 v  a  r e 


'  2  "  A  • 

These  two  values  constitute  the  two  e i genf requenc 1 ts  ; e i q e o  v a '  . e 
proper  values)  for  E  u 1 e  r i an  motion  considered  as  an  eigenvalue 
Drool em . 

For  the  first  eigenvalue  •  -  l)  ,  eq  .  ;  10-10 1  rrves 

a  =  0  and 

w  3  =  b  =  const. 

but  otherwise  arbitrary.  Thus  also 

=  .;i  -  .  =  const.  .  f  IT  9  -  N  ( 

which  gives  an  arbitra  onstant  increment  to  the  angula*' 
velocity  of  rotation  w;  ,ut  changing  the  di recti  or  of  the 
rotation  axis.  This  particular  proper  mode  r  rotation  is  the 
a  .<  i  a  1  >  p  i  n  mode  (ASM)  of  (Smith,  1  9  7  7  ). 

More  important  is  the  second  eigenvalue  -  ,  w  h 


of  :>jjrse,  i  5  tne  E  u  1  e  r  i  a  n  f  reiuency  of  polar  motion  for  one 
rigid  earth  as  discussed  in  sec.  2  .  In  fact,  new  '  10  -  10 } 
qives  a  arbitrary  and  b  -  0  ,  wnicr,  means  tnat  the  anas'' a  r 
velocity  remains  unenanged  but  the  rotation  axis  undergoes  a 
periodic  m o  t i o n  ac  cord  i  ng  to 

u  -  a  e  •  k  ,  10- 1 


.vr:  icn  is  equivalent  to  (2-16).  In  tne  terminology  of  Smite 
1977; .  this  proper  mode  is  Chandler  wobble  CW)  ,  Eulerian 
motion  being  Cnandler  motion  for  a  rig-id  earth  and  "wobble" 
b e  i  no  a  s  y  n  o  n y m  f or  polar  motion. 

Eigenvalues  for-  Spatial  Position.  The  integration  ?f 
Euler’s  equations  gives  ,  •  ,  .  ,  ., :  ,  which  define  tne  position 

O'  tne  rotation  axis  with  respect  to  a  b  o  d  y  -  f i x  e  d  coordinate 


s  /  ■  i  z  e  m 

.  Tne 

pcs  1 

t  i  o  n 

0  f 

the  rotation  ax 

i  s  wit  h 

re 

s p e c t  to 

a  n 

;  n  e  r  t  i 

a  1  s  y 

s  tern 

tha 

t  is,  its  o  r i e  n 

t  a  t  i  o  n 

a  S 

defined, 

e  .  g 

::  j  t  n  e 

Euler  angles 

:  ,  ■■ 

,.  --  requires  anotner 

i  n  te  g  r  a  t  i 

or.  , 

a,  W  e 

n  a  v  e 

s  e  e  n 

b  e  f  o 

r'  0  , 

cf  .  {  8  -  4  )  .  Let 

us  see 

now 

r+ 

-3 

(/> 

u> 

-ti 

e  c  t 

o  i  j  r  p  i 

g  e  n  v  a 

1  ues  . 

Tne 

c  u  1  e  r 

a  n  g 

I  e  s 

relate  t  li  e  bod  y 

f  ra  me 

to  tne 

i  n  - 

e  r‘  *  :  i  1 

S  7  s  t 

e  m  < 

.  f 

1 

'o  r 

tiie  present  pu*' 

nose  it 

w  i 

11  be  mor 

e 

o  n  v  p  n 

ient 

t  o  r  e 

late 

the 

body- fixed  ' 

tern  x 

X  •  X 

to  a n o 

tne 

rectangular  s/stem  xx  x  which  is  connected  to  the  inert' a", 
■/''•tern  <  i  r.  u  prescribed  simple  way.  If,  furthermore,  tne 


77 


system  x.'  is  close  to  the  body  system  x  ,  then  the  rotation 
from  one  system  to  the  other  can  be  effected  by  an  ''infinite¬ 
simal  rotation"  described  by  small  quantities  .  , 

rather  than  by  the  E  u  1  e  *’  angles  which  can  be  large. 

For  the  x/  axis  we  take  a  constant  direction 
in  inertial  space  which  is  close  to  the  figure  axis  x.  ; 
this  is  possible  since  the  x  axis  has  an  almost  constant 
direction.  Let  the  system  x :  x  _  x  .  rotate  with  respect  to 
the  inertial  system  with  cons  t  ant  angular  velocity  .'.  ,  such 
that  the  axes  x  ■_  and  x_^  never  deviate  much  from  their  uni¬ 
formly  rotating  counterparts  x  and  x 

Thus  the  system  x i  x  /  x  .  is  indeed  related  to  tne 
inertial  system  A :  X  _ X  .  in  a  prescribed  simple  way  and  can 

equally  well  be  used  as  a  reference  for  the  motion  of  the 
body  in  space. 

In  the  auxiliary  system  xi  ,  its  rotation  vector  x- 
with  respect  to  inertial  space  has,  by  its  definition,  the 
components 


iO' 
0  ' 


(  10-16) 


Since  the  deviation  of  the  frame  x 
x  is  small,  the  t ra n s f o rma t i on  from  one 
be  effected  by  a  rotation  matrix  that  is  c 


from  the  system 
to  the  other  can 
use  to  the  unit 


•n  a  t  r  i  :< 


if  the  components  of  the  v  e  c  .'  o  r  s  x  and  x  represent 
the  coordinates  of  the  same  point  in  tne  respective  systems 
The  small  matrix  ,  representing  an  "infinitesimal 
rotation",  is  s k ew - sy mme t r i c  and  may  be  expressed  as 


r  o 

7  -  -  •  .  . 

0 

Introducing  the 

vector 

(  10-13 


(  10-19  ; 


we  may  write  (10-17)  also  in  the  form 


x  =  x  •  •  a  ,  (10-20) 

where  the  cross  denotes  the  vector  product  as  usual. 

E q .  (10-20)  holds  for  any  vector  and  may  be  used  tc 

transform  vectors  from  the  x  to  the  x.  system.  In  particular 
we  have 


( 10-21 ; 


th<-'  earth's  actual  instantaneous  rotation  vector; 
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it  has,  in  the  body  frame  x  t  ,  the  components  ,  .  ,  . .  . 

entering  in  Euler's  equations  (10-2).  The  vector 
comprises  the  components  of  the  same  vector  in  the  system 
x . "  ;  it  should  therefore  not  be  confused  w  i  t  n  the  vector 

£.?  (10-16). 

In  fact,  we  have 


This  can  be  seen  in  the  following  way.  Consider-  a  point  at 
rest  in  the  body  frame  so  that  x  =  const.  then  the  differ¬ 
entiation  of  (10-20)  gives 

0  =  x_:  -  *  xj  -  ^  '  x : 

Since  the  vector  y  is  very  small  (  "  i  n  f  i  n  i  tes  i  rna  1  11  )  , 
we  shall  consistently  neglect  second  and  higher'  powers  of  it, 
retaining  only  linear  terms.  Then  the  last  term  of  (10-13)  is 
readily  seen  to  be  of  second  order  and  will  be  neglected.  Thus 
the  last  equation  becomes 

x  3  =  9  *  x 0  (10-25) 

The  comparison  of  (10-23)  with  (2-3)  shows  that  y 
is  the  angular  velocity  vector  of  the  rotation  of  the  system 
x  with  respect  to  the  system  x  ^  .  Since  the  vector  (  1 0  - 1  o  1 

describes  the  rotation  of  x  with  respect  to  the  inertial 
system,  the  sum  of  these  two  rotation  vectors  gives  the  rotation 


of  the  bo  ay  frame  witn  respect  tu  tne  merti  ai  system,  mat 
tne  earth's  actual  rotation  vector.  Tnis  proves  /  10-22). 

The  combination  of  (10-21)  and  (10-22)  now  yields 


+  1  -I'--'.  ,  (10-2 

as  usual  up  to  second-order  terms.  In  terms  of  components 
this  is 


-  h 

(10-2 

Let  us 

+■ 

now  substitute 

these  expressions 

in  Euler's 

equations  (10 

-5).  The  result 
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-  ’e)'' :  + 
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(  10-2 

=  0 


This  is  again  a  system  of  homogeneous  linear  different 
equations  with  constant  coefficients,  which  again  can  De  si  hid 
fied  by  using  complex  quantities.  Putting 


(  10-2 


+  1 


whose  roots  are 


8  2 


’  .  =  -  (10-3-lj 

The  first  two  roots  are  the  same  as  in  the  case  of  Euler's 
equations,  namely  (10-12)  and  (10-13),  defining  the  axial  spin 
mode  (ASM;  and  Chandler  wobble  (CW),  res  pec t i vel y . 

Now  the  third  eigenvalue  is  -  ,  =  -8  .  The  corresponding 
proper  mode  is  by  (10-29)  : 


w  =  te  L"t  ,  '•  .  =  0  ,  (  10-35 ) 

whereas  (10-10)  gives  for  -  =  only  the  trivial  solution 

a  =  0  -  b  ,  so  that  by  (10-7)  and  (10-9) 

-  .  =  0  . 

Tn  this  mode,  therefore,  the  rotation  axis  within  the  body 
(as  described  by  remains  unchanged,  as  well  as  the 

speed  of  rotation,  but  there  is  a  nonvanishing  w  ,  that  is, 
and  •  differ  from  zero,  which  corresponds  to  a  tilt  of 
the  rotation  axis  in  space. 

This  is  the  tilt-over  mode  (TOM)  of  (Smith,  1977).  It 
corresponds  to  a  tilt  of  the  whole  earth  (with  the  rotation 
axis  i n  v  a  r i a  b 1 y  fixed  to  it)  in  space,  so  t  n  a  t  the  earth 
rotates  with  the  same  speed  around  a  slightly  different  axis. 


8 '{ 


It  follows  that  the  TOM  Joes  not  affect  polar  'notion 
(the  axis  within  the  earth  does  not  change)  but  affects  no¬ 
tation  (the  axis  changes  periodically  in  space  according  to 
(10-35)).  It  is  thus  clear  that  the  TOM  does  not  snow  up  in 
Euler's  equation  but  appears  only  in  the  equations  (10-26) 
describing  spatial  orientation. 

Since  any  body  (regardless  of  its  internal  constitu¬ 
tion)  freely  rotating  around  a  certain  axis,  can  also  rotate 
if  the  axis  ( w ’ t  h  the  body  invariably  attached  to)  is  til  tea 
in  space,  the  TOM  must  exist  for  an  arbitrary  body,  rigid, 
elastic,  liquid,  even  inhomogeneous.  For  a  fluid  earth  model 
and  models  with  a  liquid  core,  this  mode  has  been  pointed  o u 
already  by  Poincare  '1910,  pp.  497,  50S,  513). 

These  proper  modes,  especially  CW  and  TOM,  will  play 
a  basic  role  in  the  following  sections. 


APPLICATION  TO  NUTATION  AND  POLAR  MOTION 


T n e  eigenvalue  theory  described  in  the  preceding  section 
allows  an  elegant  treatment  of  precession  and  nutation.  We 
shall  apply  metnods  outlined  by  Smith  (1977)  and  extended  by 
W  a  h  r  (19  7  9,  1930). 

We  are  using  tne  same  two  coordinate  systems  as  in 
sec.  10  :  x :  x  ;<  .  is  the  body-fixed  system  of  principal  axes 

of  inertia,  and  x:  x  -  x  .•  is  the  uniformly  rotating  auxili¬ 
ary  system.  Tne  two  systems  are  related  by 


x  =  x  + 


11-1 


up  to  terms  of  second  order  in  •  as  usual).  The  vector  _ 
describes  tne  infinitesimal  rotation  by  which  our  two  systems 
d  i  f  f  e  r  . 

The  rotation  of  the  earth  is  described  b v 


d  x. 
d  t 


al  I  vectors  refer  to  the  system  x 


3  v  ■'  10-161  and  ■'  1 0  -  i 


t  h  e  i*  o  t  a  t  i  o  n  v  e  c  tor  i 


e  :  + 


11-3 


e  -  being  the  unit  vector  uf  tin.-  .  i ...  1  . 

Since  we  ao  not  consider  changes  in  t'n;  speed  wt  r>,  t  j  t 
we  nave  „  .  -  0  •  .  {that  is,  we  0  i  s  read  **cj  toe  a  .<  i  a  :  s  p  i  n 

mode,  ASM,  see  sec.  10).  Thus 


_0 

Complex  notation  will  again  be  convenient.  We  put 

using  ■■  as  a  symbol  for  the  quantity  that  has  been  denote  a 
by  w  in  (10-27). 

We  shall  use  this  c c m p  1  e x  notation  simultaneously  wit ■' 
tnree-dimensicnal  vector  notation.  We  put 

£;=!£[,  (11-6 

that  is,  the  n  umbe  r  i  represents  rotation _ _a  round _ tne  x, 

axis  by  the  an  gle  of  +  -r  /  2  ;  e;  and  e  denote  the  unit  ve 
tors  of  the  x;J  and  x  axes.  Thus,  i  can  simply  be 

interpreted  as  a  rotation  matrix;  cf.  (Duschek  and  H  o  c  h  r  a i nor 
1961 ,  p .  222) . 

This  is  the  only  convention  needed;  everything  else 
follows  automatically.  In  particular, 


•;b 


gives  tne  connection  between  tne  comp  1  e  '•  numue 
the  c o r re s pon 0 i n c:  vector  (11-4).  in  fact,  o y 

..  -  .  0  =  ;  ■  .  +  i  •  '  e  ■ 

=  -  e  •  +  ■■  i  e  ■  =  •  •  e  +  ■  e 

identical  to  11-4).  A  relation  of  the  form  (1 
any  vector  which  has  no  comp  on ent  along  tne  x 
shall  always  use  the  same  letter  for  two  quant 
in  this  way:  the  vector  is  underlined,  tne  cor 
complex  numDer  is  not  underlined. 

Assume  now  to  be  an  exponential 

t  - -  . 

■  = 

as  in  10-29),  t  and  .  be  i  nc  constant  comp  1 e 


ce 


s  ~j  that  differentiation  is  equivalent  to  mul  fi 
i  a  fact  of  basic  usefulness,  well  known  f r  o  m  s 
:i  /  (  11-7;  , 
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si  me  relation  also  holds  for  vectors. 
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0  *  a)  Q 


-  1  a  0 


(11-20) 


and  e  is  (11-11).  Thus  (11-19)  becomes 

H  =  C(jjq  +  i  ( o  +  n )  C  e  (11-21) 

or 

H  =  Cne3  +  iA(?  +  n)  0  .  (11-22) 

This  is  in  the  body  frame.  In  the  x  ‘J  system  we  have  by 
(11-1) 

H°  =  H  +  9  x  H  .  ( 11-23) 

To  first  order,  by  (11-15), 

9_  x  H  =  1  *  Cne3  =  -iCfaO_  .  (11-24) 

From  the  last  three  equations  there  follows 

H°  =  Cne3  +  i  ( Ao  +  An  -  Cn)  9_ 

=  Cne 3  +  iA(o  -  oE)  fl  ,  (11-25) 

using  the  Euler  frequency  (10-6).  The  correspondi ng  unit  vector 
finally  is 
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•  A(o  -Or  ) 
e3  +  i  — - — 

Cii 


(11-26) 


Relation  be  tween  e_  and  torque  .  The  infinitesimal 
rotation  e_  can  easily  be  related  to  the  torque  L  of  the 
lunisolar  attraction.  The  basic  equation  is 

H°  +  *  H°  =  L 0  ,  (11-27) 

all  quantities  referring  to  the  x  0  system  rotating  with  uni¬ 
form  velocity  . 

Deriving  (11-25)  with  respect  to  time  gives 

H'  =  i  A  (  o  -  n  E  )  0_  =  -An  (  a  -  rE)^ 

using  (11-11),  and 

<  H  ’  =  ..e ->  •<  H° 

=  iA 2 (a  -  aE)e3  *  e 
=  -  i  A  ( -T  -  -E)n_  x  8, 

=  -A -  (  •  -  iE)l 

using  (11-15).  Thus  (11-27)  becomes 

'A('J  +  Q)  (  t  -  .iE)e  =  .  (  11-28) 


We  assume  L  to  have  the  form 
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L 


J 


i  lz  , 

0  i 

J 


Le  ; 


where 


(11-29) 


L  =  Lj  +  iL: 

as  before,  cf.  (3-15),  and  L  is  an  exponential 
L  =  3eiat+Y  . 

Then  the  complex  numbers  6  and  L  are  related  by 


A  (  o  +  a  )  (  a  -  a  £  ) 


This  is  equivalent  to 


e 


1 

C  12 


1 

o  +  Q 


a  -  a 


(11-30) 


(11-31) 


(11-32) 


(  11-33) 


as  is  readily  verified  by  computation. 

This  latter  form  shows  very  well  the  resonance  at  the 
proper  frequencies  -a  and  a  :  if  the  external  moment  L 

E 

has  a  frequency  o  equal  to  either  of  -the  two  proper  frequen¬ 
cies,  then  the  expression  (11-33)  will  have  a  singularity.  For 
lunisolar  effects,  whose  frequencies  are  grouped  around  the  si¬ 
dereal  frequency  Q  ,  the  relevant  resonance  is  at  nr  -  , 


! 
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which  gives  preces  s i on ;  cf.  sec.  7,  especially  eq  .  (7-14) 

(no  lunisolar  effect  would  have  a  frequency  aE)- 

Eq .  (11-32)  or  (11-33)  provide  the  possibility  to 
express  all  our  quantities  in  terms  of  the  external  torque 
L  .  For  instance,  (11-25)  and  (11-32)  give 

H°  =  Cue  3 - —  e  3  ,  (11-34) 

o  +  Q 

which  clearly  shows  the  resonance  at  a  =  -fi  .  This  shows 
the  importance  of  the  "tilt-over  mode"  (TOM,  cf.  sec.  10) 
for  precession  and  nutation. 

Nutation  and  polar  motion.  Nutation  is  the  periodic 
motion  of  any  of  the  three  axes:  angular  momentum  axis  H  , 
figure  axis  F  ,  rotation  axis  R  ,  with  respect  to  a  fixed 
reference  axis  for  which  it  is  natural  to  take  the  axis  x3° 
that  is  fixed  in  space;  it  has  the  unit  vector  e3  .  Thus 
the  nutation  vector  n_  of  any  of  these  is  obtained  from  the 
correspondi ng  unit  vector  e  by  subtracting  e3  : 

Hr  =  £r  "  (l  ’  ’ 

nF  =  eF  -  e ,  (11-35) 

(In  -  sH  -  e  , 


the  unit  vectors  being  given  by  (11-13),  (11-16),  and  (11-26). 
Using  the  complex  number  n  corresponding  to  the  vector  n 
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by 


n  =  n  e  i  , 


(11-36) 


we  thus  get 


i  -  e  > 

-is  , 


i  a(°-°e)  t3 


(11-37) 


Ca 


Polar  motion  is  treated  in  the  same  way.  It  represents 
around  the  f i gure  < 
is  defined  by  the  vectors 


motion  around  the  figure  axis  .  Therefore,  polar  motion 


p  =  e  -  e 

-R  -R  -F 


p  =  e  -  e 

t-H  -H  -F 


(11-38) 


in  analogy  to  (11-35)  (of  course,  p  =  0).  The  expressions 
(11-13),  (11-16),  and  (11-26)  give  f^r  the  co r res  pond l n g  com¬ 
plex  numbers 


.  a  +  fi  . 

p  =  !  -  e  , 


ph  =  1  Z 


A  a  +  n 


(11-39  ) 


Their  expression  in  terms  of  L  by  (11-32), 


i  L 


a“(°-3e) 


i  L 


CT (o-n£  ) 


( 11-39  '  ) 


snows  resonance  only  at  the  Chandler  (or  rather  Euler)  fre¬ 
quency  -  ,  as  it  is  natural  for  polar  motion.  The  lunisolar 

( ti dal ly-i nduced )  effect  on  polar  motion  is  not  resonant. 

These  relations  are  illustrated  in  Fig.  11.1,  in  which 
the  oriain  0  corresponds  to  the  X30  axis. 


n„ 

=  OR 

Pn  =  FR 

R 

r  R 

n 

=  OF 

F 

0  H 

=  OH 

PH  -  fh 

FIGURE  11.1  Forced  nutation  and  polar  motion  of  the 
rotation  axis  R,  the  figure  axis  F,  anc 
the  angular  momentum  axis  H  . 
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Comparison  with  previous  formulas.  In  (5-13)  we  have 


put 


L  =  (C  -  A  )  i ; :  l  B  ^  e  L  U'jt+V 


(  11-40) 


We  wri te  this  as 


L  =  v  L  . 

3  3 


(11-41) 


with 


L  =  (C  -  A  )  p. :  B  .  e  '  1  (  “jc  +  tV  •  (  11-42) 

We  shall  apply  the  preceding  results  to  each  frequency  separa¬ 
tely  and  only  at  the  end  sum  over  all  frequencies. 

Thus  we  shall  identify  L  in  (11-31)  with  an  L  as 
given  by  (11-42).  The  comparison  shows  that 

o  =  -w.  ;  (11-43) 

J 

of  course,  y  =  -  3  .  but  this  we  shall  not  need. 

Then  (11-  39')  g  i  ves 

i  L 

_ 2 _ 

A +oE  ) 

i  L. 

3 _ 

Ca(-^  +~e) 


( 11-44  ) 
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and  using  (11-42)  and  summing  over  j  we  obtain 


(  11-45  ) 

(  11-46) 


which  are  identical  to  the  forced  part  of  (5-24)  and  (6-11) 
for  a  rigid  earth  with  k  =  0. 

For  free  Eulerian  motion  we  have 


L  =  0  ,  >  =  .  (11-47) 

L 

Equations  such  as  (11-32)  or  (11-39’ )  give  0/0  and  cannot  be 
used,  but  (11-37)  and  (11-39)  remain  valid  with 

r  =  ,eiJEt  (11-48) 


for  the  proper  mode  CW  according  to  (10-29)  (f  is  the  same 
as  w).  For  free  polar  motion,  (11-39)  gives 


and  the  comparison  with  (5-24)  shows  that 

+  n 

i'T  t  E 

m  i  e  e  - , 
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so  that 


e 


O.  +  a. 


m-j  e' 


10, 


(  11-50} 


is  the  relation  between  free  polar  motion  nr  and  the  corre¬ 
sponding  infinitesimal  rotation  9  .  With  this  n  ,  the  second 

equation  of  (11-39)  gives  the  free  term  in  (6-11)  with  k  =  0. 

Nutation  is  handled  in  the  same  way.  With  (11-43),  eqs 
(11-37)  give 


n 


R»  j 


w  . 


n 


F,  j 


n 


H,  j 


-  i 


A  ( +o  ) 

3  E 


Cn 


where  ,  by  (  11-32  )  , 


9 


j 


L 

3 


A  A oj  .  (  u  +o  ) 
3  3  E 


(  11-51) 


(11-52) 


since 


a  +  n  =  -ijj  .  +  r,  =  -Aw . 

3  3 


where 


(11-53) 


as  usual. 


Substitution  of  (11-42)  and  summation  over 


gives 


n  =  i 


j  A 'Xl  .  (  *  J  ^  ) 


1J5 —  8  .  e' 1  ‘“V”6  J 


n  =  i 


i  ■>  ( ■»>  .  +  ■'  c, ) 

J  J  E 


B3e-i(V  +  V  , 


n„  -  l 


'  E 


5  4  u  .  (  \l  +  t  ~T 

J  1  E  ' 


B,e'lujttrV 


(11-54) 


To  these  forced  terms,  the  corresponding  free  terms  (11-37), 
with  -  =  ■  and  •  from  (11-50),  must  be  added.  Then  the  com¬ 
parison  with  (7-14),  (3-7),  and  (8-11)  shows  that 


i t 

+  i  '  .  s  :  n  -  -  i  n  e 


(  11-55  ) 


for  any  of  the  axes  R,  F,  or  H  ;  the  geometry  of  this  cor- 
res oondence  between  polar  motion  and  nutation  is  the  same  as 

in  ( 3  -  9 ) . 

The  factor  -i  expresses  a  rotation  by  -  9  0 L  (or  +270) 
which  has  no  deeper  significance  as  in  characterizes  only  the 
choice  of  coordinate  axes.  The  factor  e  expresses,  of 

course,  the  uniform  rotation  of  the  x;  x  j  ’  x  ;  system,  to 
which  n  refers,  with  respect  to  the  inertial  system,  to  whicn 


i 
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and  -  refer. 

In  fact,  the  present  method  offers  the  simplest  and 
most  direct  derivation  of  polar  motion  and  nutation  for  a 
rigid  earth.  For  this  particular  case  we  get,  in  a  consider 
ably  simpler  way,  the  same  results  as  with  the  approach  of 
secs.  3  through  8.  This  latter  approach,  however,  holds 
for  an  elastic  earth  and  thus  is  more  general. 

In  sec.  13  we  shall  extend  the  present  method  to 
an  earth  model  with  a  liquid  core. 
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12.  POINCARE'S  LIQUID  CORE  MODEL 


Neither  the  rigid  nor  the  elastic  earth  model  are  capable 
of  adequately  describing  earth  tides  and  nutation.  Tne  effect 
of  the  liquid  core  must  be  taken  into  account. 

The  oscillations  of  a  rotating  ellipsoidal  sneil  contain¬ 
ing  a  homogeneous  liquid  were  first  treated  simultaneously  by 
Sloudrky  and  oy  Hough  (1895).  The  most  elegant  treatment  is 
by  Poincare  (19  10).  His  paper  is  so  frequently  used  ana  :  j o  t  e  a 
that  it  iias  become  customary  to  speak  of  the  Poincare  moael. 

Since  Poincares  method  is  treatea  in  easily  accessible 
textbooks  (Lamb,  1932,  p.  724;  Melchior,  1978,  p.  122)  we 
snail  h  ere  only  describe  it  in  general  terms,  rathe  r  t  n  a  n  de¬ 
riving  it  step  by  step. 

Let  us  refer  the  ellipsoidal  shell  to  principal  a  •  e  s 
xy:  ,  then  tne  inner  ellipsoidal  surface,  whicr.  encloses  tne 
liquid- filled  cavity,  has  the  equation 


Poincare  considers  a  motion  of  the  liquid  such  that,  by  the 
t ran s f o rma t i on  (12-2),  it  is  transformed  into  a  rotation  of 
the  sphere  (12-3).  Thus,  by  (2-3),  the  velocity  in  the  auxi- 
1  i ary  x  '  y  '  z  '  sys  tern  i s 


,  x  ;■  z  *  <■>  y 

;  o  * '  -  x  i  z 

.  .  ;  y  '  -  A  c  X 


(12-4) 


if  the  co r res  pond i ng  rotation  vector  is  ^  .  Going  back  to 
the  real  system  xyz  b y (  12.2)  and  adding  the  actual  rotation 
_  ,  we  obtain 


-  5x3y  +  w;  Z  - 

•  y  . 

b 

b 

a*’x 

c  x  * L  +  J  3  x 

•JU  J  z  , 

c 

c 

pxiy 

-  gX;x  +  mY  - 

u3  j  X  . 

Here  £  represents  the  rotation  of  the  earth  with  respect 
to  the  inertial  system,  and  _<  expresses  a  rotation  of  tne  f’uid 
core  with  respect  to  the  earth.  (The  latter,  of  course,  is  a 
strict  rotation  only  after  the  formal  trans formati on  (12-2  )  to  tn 
auxiliary  x'y'z'  system,  but  for  a  nearly  spherical  earth, 
this  holds  approximately  also  in  the  actual  xyz  system.) 


The  motion  described  by  the  velocity  components  (12-5) 
represents  the  simplest  possible  motion  of  an  ideal  fluid. 

The  velocity  is  linear  in  the  coordinates  x ,  y ,  z  .  This  is 
quite  natural,  considering  that  such  a  linear  dependence 
somenow  corresponds  to  a  quadratic  potential,  and  the  luni- 
solar  potential  ii  the  usual  treatment  (sec.l)  is  indeed  qua¬ 
dratic,  namely  a  spherical  harmonic  of  second  degree. 

The  kinetic  energy  T  is  found  by  summing  (i.e.  inte¬ 
grating)  the  square  of  the  velocity  over  all  mass  elements 
d  m  : 

21  =  ///(*•  +  y  +  2 :  )  dm  . 

earth 

The  substitution  of  (12-5)  and  integration  over  the  whole 
earth  (liquid  core  plus  ellipsoidal  shell)  yields 

2  T  -  A  oj  ’  ■  +  B  -  +  C ,  j  i  ■  + 

+  A ;  -  +  B  i .  y  '  -  +  C  _  i  ,  ■  t 

+  2F^;  v  ;  +  2G... +•  2HW?X  3  .  (12-6 

Here  A,B,C  are  the  principal  moments  of  inertia  for  the 
whole  body,  and 


103 


Ac  =  }Mc(b 

-  +■ 

c’)  , 

F  = 

|-M  be  , 

Bc  =  K<c 

-  + 

)  , 

G  = 

5Mcca  ’ 

(12-7) 

C  =  iM  (a 
c  5  c 

2  + 

b2)  , 

H  = 

f  M  a  b  , 

J  c 

M  denoting  the 

C 

total  mass  of 

the 

liquid  core. 

The  equations 

of  motion 

may 

be 

made  plausible 

in  the 

following  way.  Assume 

that  the 

liquid 

core  is  absent 

and  the 

earth  is  a  simple 

rigid  body. 

Then 

the 

equations  of 

motion  may 

be  wri tten 

d  3T 
d  t  3  u  i 

3T 

i  3oj  ; 

.  3T 

+  U)  o  — - 

3aj3 

=  L. 

(12-3) 

plus  two  other  equations  resulting  from  cyclic  permutation  of 
subscripts.  In  fact,  in  this  limiting  case  we  have  A _  =  - 

=  Cc  =  F  =  G  =  H  =  0  ,  and  the  substitution  of  (12-6)  into 
(12-8)  immediately  gives  Euler's  equations  (2-12). 

It  will  now  be  assumed  that  (12-8)  also  holds  for  the 
general  case  of  a  liquid  core.  Since  the  vector  _<  plays  a 
role  analogous  to  w  ,  we  may  guess  that  the  following  equations 
also  hold: 

d  3  T  3  T  3T 

d  t  i  'X  1  (  3  3  X  2  <  2  3  x  3 


0 


(  12-9) 
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and  two  cyclically  permuted  equations.  The  change  in  sign  in 
these  two  equations  are  due  to  the  fact  that  \  ,  expressing 
rotation  of  the  core  with  respect  to  the  body,  has  a  similar 
character  as  -  ^  ,  which  describes  rotation  of  inertial  space 
with  respect  to  the  body.  On  the  right-hand  side  of  (12-9) 
there  is  zero  since  external  forces  do  not  affect  the  relative 
motion  of  the  core. 

A  rigorous  derivation  of  (12-8)  even  in  tne  presence 
of  a  liquid  core  is  not  difficult  since  it  is  simply  equiva¬ 
lent  to  the  moment  equation  (2-2);  cf.  (Lamb,  1932,  p.  724). 

On  the  other  hand,  (12-9)  may  be  derived  using  Helmholtz' 
vorticity  equation,  which  means  going  rather  deep  into  fluid 
mechanics;  cf.  (Lamb,  1932,  p.  725)  or  (Melchior,  1978,  p .  124) 

Our  plausibility  reasoning  using  arguments  of  symmetry 
between  ^  and  >_  can,  however,  be  made  rigorous.  This  has 
been  done  already  by  Poincare  (1910,  p.  484),  using  a  theorem 
on  dynamical  systems  with  groups  of  symmetry,  earlier  given 
also  by  Poincare  (1901).  This  very  elegant  theorem  is  unfortu¬ 
nately  not  found  in  standard  treatises  on  analytical  mechanics 
with  the  exception  of  (Whittaker  ,  1937,  p.  43),  (Loomis  and 
Sternberg,  196  3,  p.  541),  and  (Abraham  and  Mars  den,  1978,  sec. 
4.4),  where  similar  theorems  are  presented;  cf.  also  the  re¬ 
mark  in  (Klein  and  Sommer feld,  1910,  p.  162). 

From  (12-8)  and  (12-9),  using  (12-6),  we  derive  immedia 

te  1  y 


d  t 


f(A. 


F  ■  i  )  -  -> ;  (  B  'o ; 


Gx:  ) 


(C. 


Hv?)  =  L; 


+  G<- 


;  (  C  „  .  <-  H  .  )  f  .  f  A  J  .  f  r . 


L 

<  U  -  1  o 


r(C..  +  H  x  )  -  j  .  ( A  u .  +  f  ,  .  )  «■  ,  ;  (  B  .  +•  G  ,  )  =  L 


+  +  BcX 

d  /  r\  \  .  /It  r 


+  B..,  )  - 

•  (H.  , 

*  c  . ,  .  ) 

=■  0  , 

+  Ci;,;)  -  , 

-  ,(F. 

+  A  .  .  _  ) 

_ 

{  12-  1  1 

+  Ac--:  )  * 

,  :  (  G  „ . 

+  Bc  ..|  ) 

=  0  . 

•g-£( G-» 2  +  Bcx2  )  +  •  :  (Hu,  3  +  C 0  «.  »  )  -  ,  .  (F,  +  A  . ,  )  -■  0  , 

{  12-  1  1 

HvO  J  +  Cc  x  3  )  +■  '<  .;  (  F  .J  [  +  A  c  .  ;  )  -  ,  ;  (  G  ..  _  +  B .  2  }  =  0 

For  an  ellipsoid  of  revolution  we  have 

a  =  B  ,  A  =  B  ,  A  -  B  ,  F  .  G  (12-12 

vJ  iS 

and  furthermore,  by  (12-7), 

H  =  Cc  (12-13 

Then  the  third  equations  of  (12-10)  and  (12-11)  give 

d  /  r  .  r  \  ,  r  /  .. 


c  c  x  3 )  +  f7  ( :  X -  to  -  X  1  )  =  L  ,  , 


(  12-14) 


i 12-14 


a  . 

Jt 1 L  •  '  ' 


C  .  ,  .  )  +  F  (  ,  ; 


Xi)  =  0 


As  usual,  we  disregard  L  ,  which  causes  variation  of  rotational 
speed  Out  nut  polar  motion  (sec.  3).  Thus  we  put  Ls  =  0  and 
s  u  0  t r  a  o  t  b  o  t  n  equations  (12-14).  The  result  is  (the  dot  deno¬ 
tes  d  ‘  d  t )  : 


:ons  t .  = 


(12-15) 


nen  tne  second  equation  of  (12-14)  becomes 


C  . ;  .  +  f  : 


=  o  . 


(12-16) 


Now  ,  and  ,  have  a  similar  small  order  of  magnitude  as 

and  .  ,  so  that  ;  ,  is  a  quantity  cf  second  order, 
wnicn  w-  , nail  consistently  neglect  in  the  sequel,  as  we  did 
n  tne  preceding  sections.  To  this  accuracy,  (12-16)  reduces  to 


=  const.  ; 


12-17 


w  t  a  r  e 


i 12-13 


')  s  i  n  g  (12-13),  (12-13),  (12-15),  and  (12-13),  we  may 
write  the  first  two  equations  of  (12-10)  and  (12-11)  in  tne 
form 
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i 


A iw  ^  +■  F  ,  ■.  *  (  C  -  A  ) .  j r-  -  F  .  ^  ^  =  L  |  , 

A^  _  +  Fv  __  -  (  C  -  A )  ......  ■  +  F .. ,  ;  ■  L  ; 


(  12-19) 


F:: 

+  A  c  x  ’.  C  -  - « c 

=  0  , 

f:: 

+  Ac  +  Cc.u; 

=  0  . 

I 


Again,  complex  notation  will  be  convenient. 

u  =  +  i*.  , 

v  =  V  :  +  i  < .  , 

L  =  l\  +  iL:  , 

so  that  (12-19)  and  (12-20)  take  the  simple  form: 
Au  +  Fv  -  i(C  -  A)  .lu  +  i  F . :  v  -  L 


(  12-20  ) 


We  put 


(12-21) 


Fu  +  Acv  +  i  C c v 


=  0 


(  12-22  ) 


These  formulas  generalize  Euler's  equations  to  the  case  of  a 
liquid  core. 


Eigenvalue  problem.  We  proceed  similarly  as  in  sec.  10. 

We  put 


108 


and 


L  =  0 


Then  (12-22)  reduces  to 

A  •  i  +  F  •  ,i  -  A  t  i  +  F ti  =  0 
F  •  i  +  A.-  ’  +  Cc...  t  -  0 


(12-24 


using 


T  IT 


( 12-25 


as  usual.  The  determinant  must  be  zero : 


.A(—e)  F(:+.0  ; 

F  ■  Ac  7  +  Ccf.  ''  °  ’ 

w  n i o  h  gives 

A  (  -  -  -E)(A..  I  f  Cc: )  -  F\,(-.  +  .:)  *  0  •  (12-26 

With 

F  =  2A,.Q.  -  C„-  -  C,.  (2A,;  -  Cc)  (12-27 
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which  follows  from  (12-7),  this  becomes 

A(j  -  cE)(Ac,  +  Cc.i)  =  Cc(2Ac  -  +  ..)  •  (  12-23) 

We  now  i n  troduce 


C 

C 

the  dynamical  ell iptici ty  of  the  core.  By  (12-7)  this  is  also 
equal  to  the  core's  geometric  flattening  (apart  from  negligible 
second-order  terms ) : 


a 


(12-30) 


0(e2)  denoting  terms  of  order  as  usual. 

Then 

C  =  A  ( 1  +  e)  ,  (  12-31) 

c  c 

and 

F:  =  C  ( 2A  -C  )  =A2(l+c)(l-c) 

c  c  c  c 

=  A2(l-£2)  =a:' 

c  c 


(12-32) 


1  10 


since  terms  of  second  order  in  will  now  consistently  oe 

negl ected . 

Thus  (12-28)  becomes 
AAc('  -  .:e)(-j  +  +  ■--)  =  Ac-a(;  +  :.) 

o  r 

A.;  -  (C  -  A):; :  (  i  +  +  £;  )  =  A  i  )  .  (12-33) 

C 

The  introduction  of  the  principal  moments  of  inertia  for  the 
mantle 

A  =  A  -  A  C  =  C  -  C  ,  (12-34; 

me  me 

allows  a  reduction  to  the  form 

(A  .i  -  (C  -  A ):.:](.-  +  il)  =  -[Ao  -  (C  -  A ) ..  ( .1 c  .  (12-35,. 

rn 

This  equation  is  solved  by  successive  approximations. 

We  first  put  =  0  .  Then  the  right-hand  side  is  zero,  and 

the  equation  has  the  two  roots 

j  C  -  A 

«.T  ,  =  -  ,-i 

,  at  > 

A  m  (12-36) 


As  a  second  approximation  we  put 


For  j  =  ,3 ;  we  have 


A  ?, 

m  * 


( C  -  A ) ;:  =  A  ci  1  -t-  A  r  i  -  (C  -  A).. 
'  m  •  m  ' 


so  that  ( 12-35  )  becomes 


(12-37.) 


£  -  i  Am(o  i  "  +  2)  =  -  A ( ■,:  t  -  a£)w.;  . 

Now  e  is  the  core  flattening  (12-29).  Both  ,  by  (  12-36 

and  ,  by  (12-25),  have  the  same  order  of  magnitude  as  ■.  . 

Thus  the  right-hand  side  has  actually  the  order  of  c  and 
can  be  neglected.  Hence,  to  our  usual  linear  approximation 
we  have  t  =  0  and 


T.  (  12-3S) 

n 

This  is  the  frequency  of  the  Chandler  wobble  (CW)  for  the 
present  case.  It  differs  from  CW  for  the  rigid  body,  (10-6), 
by  the  replacement  of  A  in  the  denominator  by  A  ,  the 
principal  moment  of  inertia  for  the  mantle. 

A  new  feature  is  brought  into  the  picture  by  the  first 
root.  For 


we  have 


r 


so  tnat  (12-35)  becomes 

A.n  ?  -  (C  -  A  ).)■::  ■  -  -  .  An ,  :  -  (C  -  A) 

or  with  (12-36), 

-  A,,.  +  C  -  A )  :  =  ( A  +  C  -  A ) ..  , 


w  i  r.  n  the  solution 


C 


A  r  C  -  A 
m 


( 12-39) 


Now,  using  (12-31)  and  (12-34) 

A  +  C-  A  =  C-  A  =  C  -  C  (  1  -  e  ) 

in  c  c  '  1 

=  C  -  C  +  0  ( £  )  =  C  +  0  ( e  )  . 

Since  is  multiplied  by  c  in  (12-37),  the  term  0  ( ) 

is  multiplied  with  >;  to  give  a  negligible  second-order  term. 
Thus  (12-39)  becomes 

C 

r-  =  "  C 


( 12-40) 


1 1 J 


and 


in 


This  root  has  no  equivalent  in  the  rigid-body  rotation. 
In  particular,  it  is  not  the  direct  equivalent  of  the  tilt  - 
over  mode  although  it  is  numerically  very  close.  We  snail 
come  back  to  this  question  at  the  end  of  the  present  section. 

To  give  an  idea  of  the  order  of  magnitude,  we  take 
(Melchior,  1978,  pp.  129) 


A  C 

s  0.11  =  -£■ 

A  C 

.  1 

*■  '  ttotj  • 


(  12-42  ) 


Then 


C  l  i  1  1  o 

IT  '  T'"-  C  /C  ~ 

m  c 

‘f  '  tot  *  °-0028  • 

m 

so  that 


■j-,  =  -  1.0028 a  . 


(  12-43  ) 


i 

1 


This  is  tne  eigenvalue  for  t n e  n early  cl _u r n a  1  free 
wobble  (NDFW)  according  to  the  terminology  of  (Smith,  1977, 
would  represent  a  diurnal  period). 

Tne  comparison  of  the  present  treatment  with  sec.  10 
shows  one  comparable  (though  not  numerically  equal)  eigen¬ 
value,  namely  the  C  W  value  -  •  .  Our  present  _  ,  for  the 
NDFW  ,  has  no  equivalent  in  the  case  of  ruler's  equations; 
an  tne  otner  nand,  we  have  not  yet  obtained  the  eigenvalue 
C  corresponding  to  tne  axial  spin  mode  (ASM)  .  Its  pnysica! 
interpretation  --  rotation  about  the  same  axis  witn  a 
slightly  different  speed  of  rotation  --  shows  that  it  must 
be  quite  general  and  not  restricted  to  a  solid  body.  In  fact 
we  have  it  also  in  tne  Poincare  model.  The  equation  (12-15) 

=  0 

nas  the  exponential  solution 


w  i  t  h 

i  ■  .  e  ’  -  0  , 

which  aives 

(12- 


0 


'  or 


as  another  eigen v a 
in  sec.  10.  In  the 
also  from  '12-17). 


u  e  .  This  A j  M  v  a  1  u  <  • 


no to 


present  problem,  it  >  j  u  1  d  be  o  b  t  a  1  n  e  a 


Spatial  position  .  S  o  f  a  r  ,  our-  e  .<  p  o  s  i  1 1  o  n  nas  c  j  r  - 
responded  to  the  treatment  of  Euler's  equation  in  sec.  10, 
so  that  we  have  not  yet  obtained  the  tilt-over  mode  T  .1  ’-1 
which  should,  however,  occur  in  the  present  model  as  we1!, 
for  similar  physical  reasons  as  the  ASM  . 

We  therefore  proceed  as  in  sec.  10,  considering 
spatial  position.  Since  the  mantle  is  rigid,  the  reasoning 
leading  from  (10-17)  to  (10-2  5)  holds  uni. hanged  for  the 
Poincare  model,  and  the  first  two  equations  of  10-25;  can 
be  combined  in  complex  notation  to  give 


u  =  w  +  i w  , 


i  2  -  4  0 


using  (10-27)  and  (12-21).  The  substitution  of  tin's  into 
(12-22),  and  of  the  third  equation  of  (10-25)  into  (12-15 
gives 


A  ( w 

+  i  f.w ) 

-  i  (  C  -  A  )  w  <  w  +  i  ,.w  ) 

F  ( w 

+  i  ~i  w  ) 

+  Av+iC.:v  =  0  , 

c  c 

•3  > 

=  0  . 

+  F  (  v 


..v)  -  L  , 

(  12-46 


The  exponential  form  (proportional  to  e '  t)  for  v,  and 
3.  leads,  as  usual,  to  a  system  of  homogeneous  linear  equa¬ 
tions  (if  we  put  L  =  0).  The  condition  of  vanishing  deter mi- 


i  ;  o 


nant  may  be  brought  to  the  form 

A  (  -  -  ■ ...  )  ( A  .  •  t  C  ...  )  -  F  •  (  '  +  '•).(■  +  J  =0  . 

'  12-47; 


The  solutions  of  this  equation  are 

•  -  =  0  (ASM)  , 

C  -  A 
.  .  =  — 

•.  =  --(I 

For  •  and  •  this  is  immediately  obvious,  and  and 

are  seen  to  be  the  same  as  (12-38)  and  (12-41)  by  noting 
that  the  expression  between  parentheses  is  identical  to  (12-26) 
Although  the  modes  TOM  and  NDFW  are  conceptually 
completely  different,  their  numerical  closeness  is  so  striking 
that  it  is  tempting  to  look  for-  a  physical  interpretation  of 
N 0  F W  in  terms  of  TOM.  A  hint  is  provided  by  noting  t  n  a  t  for 
a  strictly  spherical  core,  with  -  =  0  ,  we  have  :  :  =  -  . 

Thus,  for  =  0  ,  NDFW  coincides  with  TOM. 

Now  the  tilt-over  mode  characterizes  a  tilt  of  the  body 
with  respect  to  some  external  reference.  In  TCM  in  the 
proper  sense,  such  an  external  reference  is  inertial  space.  An 
ideally  fluid  spherical  core  is  mechanically  completely  inde¬ 
pendent  from  the  mantle  (''decoupled")  since,  in  the  case  of 


(CW)  , 
(NDFW)  , 
(TOM)  . 


f  12-48 


it 


m 


spherical  symmetry,  a  coupling  could  only  be  effected  by 
drag  of  friction,  which  is  absent  with  an  ideal  fluid.  Thus 
the  core,  being  independent  of  the  mantle,  can  serve  as  an 
external  reference  for  TOM  into  which  thus  NDFW  degenerates 
for  a  spherical  core. 

This  decoupling  is  no  longer  true  if  the  core  is  ellip¬ 
tical.  Due  to  the  unsymmetry,  there  is  now  a  mechanical 
coupling  :  the  inertia  of  the  core  resists  a  rotation  of  the 
shell.  Thus  we  have  an  i nertial  coupl i ng ,  or  Poincare  coupling, 
between  core  and  shell  which  is  zero  for  •.  =  0  and  can  be 
expected  to  be  proportional  to  t-  by  small  £  .  This  is  indeed 

borne  out  by  (12-48)  :  the  deviation  of  .1  s  from  is 

proportional  to  e  . 

A  detailed  study  of  the  mechanical  situation  from  a 
somewhat  different  angle  is  found  in  (Toomre,  1974). 

Proper  modes  are  also  called  resonant.  The  presence 
of  two  different  but  almost  equal  eigenvalues  0 :  and  - 
causes  a  significant  deviation  of  the  rotational  behavior  of 
an  earth  with  a  liquid  core  from  that  of  a  rigid  body,  which 
by  Poincare  (1910)  has  been  called  doubl e  reson a  n  c  e  . 


1  1  Q 
i  lu 


13.  LIQUID  CORE  EFFECTS  ON  POLAR  MOTION  AND  NUTATION 


Let  us  first  introduce  a  convenient  termi no  1 ogy  : 

u  =  jo  ■_  +  i ...  _  .  body-referred  rotation, 

v  =  >;  <•  i  \  '  .  core  rotation, 

w  -  ;  *  i-.  =  v  .  space  rotation. 

In  fact,  ( _  ,  .,  )  express  the  position  of  tne  actual  rotation 

axis  with  respect  to  tne  body-fixed  system  of  figure  axes; 

(  ;  ,  .-)  are  a  measure  of  the  rotation  of  the  core  with  re¬ 
spect  to  this  body  frame;  and  (■;  ,  ,  )  characterize  the  rota¬ 

tion  of  the  body  with  respect  to  inertial  space.  The  quantity 
;  +  i-  has  been  denoted  by  in  secs.  10  and  11  and  by 

w  in  sec.  12;  we  shall  continue  to  use  this  notation. 

The  body- referred  rotation  u  characterizes  polar  motion. 
In  fact,  the  complex  number  m  describing  polar  motion  is 
related  to  u  by 

u  -  m  ;  (  13- 1 ) 

c  f .  (3-13)  and  ( 3  -  15  )  . 

The  quantity  u  is  a  solution  of  the  basic  equation 
(12-22).  The  eigenvalues  of  this  equation  have  been  found  to 
be  ( 12-38) , 

C  -  A 
'  ~  A  " 

m 


(13-2) 


for  CW ,  and 


-  =  ‘(I*  f  )•:  (  13-3; 

-  m 

for  NDFW. 

Thus  the  general  free  solution  (without  external  forces) 
has  the  form 


The  first  term  constitutes  the  usual  Chandler  wobb 1 e ,  the 
second  term  is  the  nearly  diurnal  free  wobble.  The  principal 
contribution  to  polar  motion,  of  course,  comes  from  CW,  and 
it  is  an  open  question  whether,  for  the  real  earth,  the  NDFW 
has  a  coefficient  which  is  significantly  large  to  be 

observable  at  present  (Rochester  et  a  1 . ,  1974  ;  Yatskiv,  1980). 

Practically  much  more  important  is  the  effect  of  the 
NDFW  on  forced  motion.  Consider  the  inhomogeneous  equation 
(12-22),  the  torque  L  being  given  by  the  exponential  (11-31), 

L  =  8ei(ot+Y)  ,  (13-5) 

which  represents  a  typical  term  in  an  expansion  such  as  (5-12) 
with  a  =  -j  .  .  We  put 


u  =  Ue 


v  =  Ve 


'  15-7 


The  substitution  into  (12-22)  with  F  =  A  by  (12-32),  tnen 
leads  to 


i  A  ’  U  +  i  A  _  •  V  -  i  (  C  -  A  )  .  U  +  i  A  ..  V  =  i  K  , 
i  A  '  U  +  i  A  -  V  +  i  C  ..  V  =0 


A,-  -  -  A).:;  U  +  A  )  V  =  K  , 


A  y;!J  f  (A  1  +  C  ;; )  V 


=  0  . 


'he  solution  by  means  of  determinants  gives 


A  -  +  C  ,T 


(13-9;, 


( 13-  io ; 


where 


is  the  determinant 


A'<  -  (C  -  A ) )  ( A  -7  +  C  T )  -  A  •  (a  +  y )  .  (13-11) 


Of  principal  interest  is  U  which  gives  polar  motion. 
Let  us  compare  U  with  the  value  which  would  correspond 
to  the  same  moment  K  if  the  earth  were  rigid.  Then  Euler's 


equations  or  simply  (13-8)  for  A  -  0  C  give 


Ao  -  (C-A).: 

From  (13-9),  (13-11),  and  (13-12)  we  thus  get 


..  Ao  -  (C-A)3  i  (A  o  +  C  ..)  -  A  ■  ■  ■  (  ■  +  .) 

Uq  _ _ c  c _ 

^  [  A  ,7  -  (C-A),;  j  (A  ’  +  C..,0 

A  ^;(a  +  .1) 

=  1 - - -  .  (13-13) 

[Ao  -  (  C  -  A  )  ft  ]  (A  o  +  C  .0 
c  c 

This  can  also  be  brought  into  the  form 


U>  - 
U  ' 


1 


A  o(o  +•  ft) 

_ c _ 

A  (  0  -  n  E  )  (  o  +  ft  +  e  ft  ) 


(  13-14) 


where  o  is  (12-25)  as  usual. 

If  we  change  over  to  our  usual  notation  for  polar 
motion,  precession,  and  nutation,  we  must  put 


o 


*  'Jj 


j 


•  A  fjj  .  _  Q  * 
D 


A  uj  .  =  w  .  -  f.  =  -  ( o  +  ft ) 


(13-15) 


Then  (13-14),  on  putting  q 


U  / U  ,  becomes 


||  “  1 

q  =  -  =  ,  1  - - -_J - i - ■  ,  1  13-16 

Un 

-  A(.?  ♦  ,E)(^J  *  -)  • 

identical  to  eq.  (6-43)  of  (Melchior,  1978,  p.  128); 
note  that  Melchior's  .  has  opposite  sign. 

Thus  the  amplitudes  of  forced  polar  motion,  computed 
for  a  rigid  earth,  must  be  multiplied  by  this  factor  to  get 
the  co r re s pond i nq  amplitudes  for  Poincares  model.  Later  in 
this  section  we  shall  see  that  the  factor  U  /  U  ->  also  holds 
for  the  amplitudes  of  nutation  of  rotation  and  figure  axes 
(the  nutation  of  the  angular  momentum  axis  is  the  same  as  for 
a  rigid  earth). 

Numerical  values  are  seen  from  the  following  table 
which  is  taken  from  (Melchior,  1978,  p.  129). 


Precession 

Principal  nutation 

(18.7  years) 

Annual  nutation 


i  q  =  U^/U 


0 

1 

e  u  /  6  8  0  0 

0.994 

-m/6800  : 

1.007 

+  >i  /  3  6  5 

3.458 

-m/365  1 

1.062 

- m/  183  | 

1.260 

Semiannual  nutation 


Fortnight’ y  nutation  * . ;  1  j .  /  1  .  1  '  1 

-  .  /  H  .  7 


This  facto'*  thus  remains  1  for  p  re  l  e  s  s  i  r.  n  ,  out 
becomes  >  for  the  proper  f  req  uen  c  i  o-.  and  3"  * 

oe  s  n  ow  n  later-. 

Angular  momentum  and  to  ''sue  .  T  n treatment  f 
is  easily  extended  to  the  Po  i  ncare  model  uiu-  the  mj  n  t  i  a 


still  rigid  in  this  model  a  no  the  mode  ’  i  axially  symmetri 
Therefore  tne  Kinematics,  up  to  11-  16),  rein* i ns  uncharged . 
What  changes,  is  the  dynamic:,. 


The 

angular 

momentum  e 

X  p  r  e  s 

u  i  o  n  , 

1  1  - 

1  /  i  is  'low  re ',)  1 

by 

H  = 

C  ^  +  D 

1  } 

where  C 

is  given 

by  (  1  1  -  1 

a  n  d 

D  ’  s 

the 

m  a  t  r  i  x 

'  r  C 

0 

D  = 

0  G 

0 

;  1  3 

0  0 

H 

Tne 

quanti 

ty  D 

is  an  a d d i r 

i  o  n  a  1 

annul 

a  r 

momentum  due  to 

the 

motion 

of  the  ' 

liquid  rore 

;  L  .i  mb 

1  Q  3  3 

,  P 

.  7  34  ,  eq  . 

it  is  indeed  a  unit  vector  since  e  .lifters  from  1  on 


by  small  terms  of  second  order  in 

To  find  tne  relation  between 


a  n  d  v  . 

and  tne  annular  m o 


mentum  j_  ,  we  may  proceed  as  follows.  Equation  12-4 
exponentials 


u  -  U  e 


w  -  We 


becomes 


U’6 


T  h  u  s  j  a  /  >  i  i  -  3  2  ) , 

L 

w  -  r,  w  =  -  q  — - — - - — — - -  .  ,13- 

A(  ■  +  ■■){'  -  1  £  ) 

Inis  is  the  desired  relation  between  and  L  .  It  snows 

t  n  a  t  tne  factor  q  ,  eq.  (13-16)  which  converts  a  from 
the  r'lqid  case  to  the  Poincare  model,  likewise  converts  w 
It  could  appear  that  w  becomes  singular  if  -  =  •_ 
although  h n i s  is  not  an  eigenvalue.  Such  a  singularity  is 
nut  real,  however,  since  U  =  0  for  -  -  •  .  This  becomes 

lean  ;f  we  e<pre;s  u  by  (13-9)  ;  by  the  a  i  a  of  ,13-6,. 

\  n  j  1  •  -  2  3  t  n  i  s  becomes 

A  •  +  C 

,i  -  -i  - - —  L 

'  i »-  “  e  '"iii  t.  a  r.  t  •  can  obviously  be  written  in  the  form 


•  >  r  •  n  ••  r  i  he  ri  v a  i  ue  s  •  and 

ii  . . 


and  (13-2-1)  gives 


A  A  ( 

m  c' 


)( -  -  \  ) 


An  expansion  into  p a  r  1 1  a  1  fractions  gives  expressions  j f  t n o 
form 


{  1  j  -  J  u 


jinere  ^  and  2.  are  constants  arid  the  eigenvalues 
• ;  are  given  by  (12-48).  Thus,  polar  motion  u  is  resonant 
for  CW  and  N D F W  ,  and  space  rotation  w  is  resonant  for 
C  W  ,  N D F W  ,  and  TOM  ,  as  it  should  be. 

The  relation  between  H  and  L  is  found  as  follows. 

W  i  t  h 

_f  =  •)£:  =  we;  =  -  ■  o  •  ,  (  13-31 

by  ,'13-24',  we  may  transform  (12-21'  as  follows 


r >  e  first  i1 1  j  u  a  r.  i  o  r.  of  ;  1  3  -  s  ,  ,  o  n  n  i  t  i  o  1  y  i  n  g  ay  s  * 
i  'i t  >  r, !  n  .]  i  a  r.  o  a  c  c  a  unt  (13-6  }  ana  ;  1  3  -  7  ;  .  ns  a  y  be  w  ’’  i  ttei 

A ;  -  • o  +  A  (  -  a  )  v  -  -  l  L 


mb’ ring  2 1 -.?■<? )  and  '  13-33)  ana  noting 


,  we  get 


i  i  n  is  tne  ••  »me  relation  as  for  a  rigid  oody.  Tni  s  con  - 

r"i  .  t  re  i  ride  pen  don  e  of  tne  '-elation  between  angular  mo- 

r  and  t  > r • ; u *?  L  from  the  interna:  structure  of 
:  o  ::  o  d  ,  ;  :f‘.  t  n  e  •:  0  n  c  !  u  d  i  n  q  remark  of  sec  .  . 

hit  at  ion  and  po i a  r  motion.  Now  it  is  s t ra i gn t f o rwa 

■■  tne  "  o  rm:j  1  a  s  for  precession  '11-37;  and  polar  mo 

*  a  tne  present  case.  Since  the  first  two  relations 
■  >:■  p tire  1  y  i  n e m a  t  i  c a  1  ,  tney  continue  to  hold  : 


n 

'  -  l .  -  1  ,‘J 

I  V  °  1 

,  \  h  o  n  i  i  t  a  t  i 

o  n 

1  1 

K  .  7  h  ^ 

P  re 

::  wj  . 

vain-1  w 

o  V 

w  ,  of  c: o u r s e  ,  differs  from  tne 

factor  q  ,  according  to  '13-35 


Thus  the  amp  1 i t  u  J  e  s  of  n  u_t  a  ti  o  n  to  r _ t  h  e  ro  t  a  t  i  o  n  axis  m  d  r.  n 

figure  axis  d  i  t  e  r _ from  the_  correspond  1  ng  rigid-body  /a  lues 

by  the  factor  u  . 

This  is  not  true  for  the  nutation  of  the  angular  niomen 
turn  axis.  From  (13-22)  and  (11-3 o)  we  get 

A(:  -  :£)  A.. 

n ..  =  i  - —  w  +  v  i  12-3 

H  r  -  r 


which  differs  from  (11-37)  by  a  term  due  to  core  motion. 
Expressed  in  terms  of  l  ,  we  even  gcj.  exactly  the  same 
nutation  as  in  the  r i g i d  bod yeas e  : 


(  12-3 


by  (13-34). 

The  polar  motion  of  the  rotation  axis,  being  purely 
kinematical,  remains  the  same  as  (11-39)  : 


.  +  d  u  ,  ,  ,  , 

p  =  i : w  =  -  (  13-a 

R  i i 

and  consequently  differs  from  the  rigid-body  case  by  the 
factor  q  . 

The  polar  motion  of  the  angular  momentum  axis  is  dif¬ 
ferent  from  (11-39)  : 


V 


f  13-40} 


P., 


Here  y  obtained  from  (13-10),  which  by  (13-6)  and  (13-7) 
becomes : 

A  - 

v  -  i  — —  L  (13-41) 


These  relations  are  illustrated  by  Fig.  13.1,  which 
is  the  extension  of  Fig.  11.1  to  the  Poincare  model.  It 
shows  that  the  rotation  axis  R  is  no  longer  close  to  the 
angular  momentum  axis  H  . 

Free  motion.  In  the  absence  of  external  forces,  for 
L  -  0  ,  the  basic  equations  (12-22)  have  the  solution: 


A  .  t,  +  (  A  .  +■  C  ‘-J )  2.  .  -  0 

Oil  '  1  c  1 


:  +  (1  + 


.  +  (1  +  -J: 


w  i  t  h 


E  r;  s  .  (13-44)  are  satisfied  if 


;)'■  , 


1  =  -  <  ’  ' 


13-44 


13-45 


13-46 


with  arbitrary  complex  constants  *  and  u  ;  this  is  easily 
seen  to  be  the  general  solution.  Thus 


u  =  (1  +  -  ;  ) '  e  " 1  +  (1  +  -  .  )  u  e : 


v  -  -  •  i  '•  e 


(  13-47) 


represents  the  general  free  solution. 

Important  is  the  order  of  magnitude.  From  (13-45)  with 
(13-3)  and  (12-43)  we  get 


1  3  3 


C  -  A  =  C  -  A  .  C  -  A  .  1_  / 

A  A  -  A  A  300 

m 


and  similarly  using  (13-3)  and  (12-31) 


1  a  ^ 
1  +  £  Q 


1  +  »-c/cm  _ 

1  +  t 


=  -  (  1  +  ^  -  «:  ) 
'■'m 


c, 

c 

m 


(  13-48) 


Hence,  approximately, 


u  =  \ei0lt  +  0  (  c  )  u  e  1  a  •'  t  , 
v  =  O(e)>ei0lt  +  ye  10 2 1 


(  13-49) 


This  shows  that  polar  motion,  (u)  comes  principally  from 
(frequency  a;  )  but  core  motion  (v)  is  caused  mainly  by 
( frequency  o2  )  . 

The  corresponding  free  value  of  spatial  rotation 
is  given  by  (13-24)  where  a  =  o  \  or  ,  depending  on 

mode.  Thus  from  (13-47)  we  get 


CW 

NDFW 


the 


134 


using 
in  the 

Tnen  , 

and  us 


We  re c 


13-50 


Now  free  nutation  can  be  computed  by  (13-35)  and  (13-36) 
(13-4?)  and  (13-50),  and  free  polar  motion  is  obtained 
same  way  from  (13-33)  and  (13-40). 

It  is  very  instructive  to  consider  the  N D F W  only, 
putting 


re1"--  =  h  (13-51) 

ing  (13-48)  and  (13-3)  we  get  from  (13-47)  and  (13-50)  : 

C.  . 

u  ------  n  , 

c« 

v  (1  +  --)h  ,  (13-52, 

"  i:i 

,  cc  . 

i  w  =  —  n 

C  . 

all  that  the  core  flattening 


1_ 
4  00 


and  that  for  the  me  merits  of  inertia  of  core  and  mantle  we 
h  a  v  e 


1 


c  c  c  -  c 

—  =  0.1  ,  —  =  - —  -  ■-  0.9 

c  c  c 


This  shows  that  the  ratios 


m 

C  1 
-  _jh 

C  c 

c 

i  ,lw  C  1 

___  -  -  —51 

u  C  e 

are  both  very  large. 

In  particular  we  have  for  this  mode  by 

(  13-40) 


p  C  C  9  p  C  Cu 

R  R 


A  A  C  1 

•  c  m 

C  c  C  C 


(13-38) 


The  first  term  on  the  r i g h  t - h a n o  s . d e  i  ;  close  to  :  a n o 
can  be  neglected  with  respect  to  the  very  large  second 
term.  Furthermore,  A  /C  =  1  r  0 ( ►  )  .  Thus 

P  C 

_J±  =  _  _i£  =  -  400  .  (  1  3-53  ) 

Pr  C, 

Using  (13-3)  this  can  also  be  written 


This  simple  relation  permits  a  kinematics!  i n te r p re ta t i on  in 
terms  of  body-  and  space-fixed  cones  according  to  Poinsot 
(Rochester  et  a  1  .  ,  1974). 

The  nearly  diurnal  free  wobble  is  illustrated  by  Fig. 
13.2,  which  in  view  of  (13-53)  shows  that  the  angle  a  between 
the  angular  momentum  vector  and  the  figure  axis  is  about  400 
times  larger  than  the  angle  -  between  the  rotation  axis 
(vector  .>)  and  the  figure  axis.  Cf.  also  (Toomre,  1974  ;, 
(Rochester  et  a  1  .  ,  1974)  and  (Yatskiv,  1930);  Rocnester  et  a i  . 


FIGURE  13.2  The  nearly-diurnal  free  wobble. 

which  differs  from  (13-53)  only  by  0  (  .  )  which  we  have  disre¬ 
garded  . 

The  same  phenomenon  may  also  be  looked  at  from  a  si  ignt- 
ly  different  angle.  In  the  absence  of  external  forces,  the 
angular  momentum  H  retains  its  position  in  space  j  the  nuta¬ 
tion  of  the  angular  momentum  axis  is  zero  : 

n  =  0  .  (  13-56) 

n 

This  is  evident  from  (13-37)  since  the  denominator  differs  from 
zero  if  o  -  •  or  a.  ,  and  the  numerator  is  zero  if  L  =  0  . 


,  >  '  ”  1  *  i  ,  • 

n  ' 

-f '  ■  4  :j  i  j  .  1  ;  - :  ■: 

p  _ 

T  o  f.  he  near  i  y  diurnal  r  r  e  <•  w  •  >  r>  b  1  e  there  exists  a  comet  pone  i  n  c 

nutation  which  is  ahou  t  4  0  0  times  larger.  Such  a  nutation  na 

not  yet  been  observed,  wnr.ti  is  anotner  indication  that  tie 
a  mu  I  i  t  tide  of  tne  MDFW  for  tne  eartn  must  be  very  small. 

Compilation  of  formulas.  The  basic  quantities  for  f re 
motion  are  body  -  re f e r red  rotation  u  ,  core  rotation  v  ,  and 
space  rotation  w  ,  which  we  collect  first: 

i  _  -  (1  *  •  in'  •  '  *  ;  1  +■  •  .  )  ..e  ‘ 

v  - •  .  • e '  -  ■  ■  r e  ,  i  13-59 


Then  the  tree  polar  motion  of  the  rotation  axis  R  and  the 
angular  at  omen  t  um  avis  li  become 


'lo'Vtr  n  t  uin  a 


-  :  r_  *1  e  a r1  a r  a  r 


Forced  m  o  1 1  n  *'o  r  th*f  a  n  <5'.j  i  .a  r 
same  a  ‘  for  the  rigid-body  case ;  fo  r  tn 
factor'  o  as  i  13-35)  shows.  I  h  u  s  t  m? 
sec.  11,  e  o  s  .  1  1  1  -  5  a  )  and  (11-50),  g i  v  e 

-  i  .  .  •  f.  •  i 

;  1 

^  i  .  s  ’  n  )  r  '■  '■  - 


Toe  total  polar  mot'ur-s  m«i  nut  a 
a-j  :  >  r  3  *  o •.*  c  r  re  s  o o  o d  :  n  o  1  »'<■*.;  :  ml  t  o  rc 

o  r 0  V ' ded  h /  t h e  r a ; t  foot,  (or  -  ■ . 

an.j  q  -  1  ,  nest  dive  the  .  a.ne  riqid- 

5  _  2  ,  .o-ll',  •  -  1  J  1  ,  /  :~i-  '  i  ,  .»  nd  •  ••;-  • 

precessions!  par-"  pro  nor  t  ’  wa  !  to  LI 

•'em  airs  t  ?i  s  a  me  a  1  s  ...  t  o  >■  i  n  e  p  r  o  s  e  n  t 


s  r  ->  <>  i  r. •:  t.  i  on  o-i '.a  i  *  n ••  forced  nutation  of  Vie 

•  I  ••'-O'  ’  ;  -  i  ♦  1  '  .COS’’)'  ’  ‘  .  '  13-  ' 

/  p 

i''  1  u '!  i  n  renia  rf  .  An  adequate  e  a  r  t  '■  mode'  must  ta*e 
ounf.  ■  m?  e  1  i '■  ticify  of  tne  mantle,  3  j  well  as  a  core 
o 1 1  f  >n  ileal  fluid.  1  h  is  co  mp  1  e  .<  and  difficult  t  c  o  1 
t  .  1  o  i  rp.  .i  re  i  19  10  ■  and  nas  oeen  extensive1)'  stjdie 
1  i  -  f  decade1'  ,  hi  anil  niits  peine  •'  .Je  -  -  reys  and  V  i  c  e 
'■!  .  a e n  •.  r  ,  10  6  1  i  ,  i  nd  W a  n  r  ,  1 9  '  9  •  .  Tne  c  e  ne  t  n  o  a  s 

of  I  *  *  i.  t  ci  i  a  no  t  h  e  r  re ;;  v  ••  t . 
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